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INTRODUCTION 

The  feature  of  mathematics  which  we  may  call  scientific  approxi- 
mation did  not  receive  the  recognition  among  the  ancient  Greek 
geometers  that  has  been  accorded  it  by  mathematicians  of  modern 
times.  Since  their  work  in  pure  mathematics  was  based  on  geometry, 
the  Greek  scholars  evaluated  irrational  roots  as  readily  as  rational 
roots  by  the  construction  of  lines.  But  when  it  came  to  applied 
mathematics,  numerical  approximation  was  often  not  only  desirable 
but  absolutely  necessary;  and  so  we  have  some  remarkably  close 
approximations  by  Archimedes  (c.  225  B.  c.)  in  his  work  in  mechan- 
ics, by  Heron  of  Alexandria  {c.  200  A.  D.)  in  his  work  as  surveyor, 
and  by  Theon  of  Alexandria  (c.  375  a.  d.)  in  his  astronomical  com- 
putations. Archimedes  even  used  numerical  approximations  in  the 
domain  of  pure  geometry  in  his  Measurement  of  the  Circle.  On  the 
whole,  however,  the  Greek  mind  did  not  take  to  this  form  of  mathe- 
matical activity. 

The  Hindus  and  the  Arabs  used  approximation  methods  more 
extensively  than  the  Greeks,  and  notable  work  was  done  by  Brah- 
magupta  (c.  510  a.  d.),  Bhaskara  (c.  1150  a.  d.),  and  an  unknown 
Arab  computer  mentioned  by  Chelebi.  However,  it  was  during  the 
Renaissance,  when  the  newly  discovered  general  formula  for  solv- 
ing the  cubic  was  found  to  be  inoperative  for  the  "irreducible  case," 
that  approximation  methods  became  a  vital  problem,  and  for  two 
centuries  the  keenest  mathematical  minds  worked  to  find  smooth, 
effective  methods  for  approximating  cube  and  higher  roots. 

The  purpose  of  this  research  is  to  trace  the  history  of  the  different 
methods  of  approximating  roots  of  numerical  higher  equations  that 
were  used  up  to  1819,  the  date  of  the  publication  of  Horner's  method ; 
we  purpose  to  trace  their  early  beginnings  in  finding  the  roots  of 
numbers  and  in  solving  incomplete  equations,  and  to  watch  their 
growth  into  systematically  developed  general  methods  for  solving 
complete  equations  of  any  degree.  Pure  trial  and  error  methods, 
like  some  of  the  methods  of  the  Egyptians  and,  in  modern  times, 
that  of  Junge,  will  not  be  taken  up,  nor  shall  we  discuss  such 
ephemeral  methods  as  those  of  Fontaine,   Collins,  and  some  of 
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DeLagny's.  We  shall,  moreover,  limit  the  investigation  to  algebraic 
methods  in  algebraic  equations  of  one  unknown,  and  only  inci- 
dentally shall  we  touch  upon  transcendental  equations  and  geo- 
metric and  trigonometic  methods. 

In  searching  for  literature  on  this  subject  the  author  found  only 
two  articles,  outside  of  the  encyclopedias,  in  the  nature  of  surveys: 
one  by  Augustus  De  Morgan  in  the  Companion  to  the  British  Almanac 
(1839),  entitled  "Notices  of  the  Progress  of  the  Problem  of  Evolu- 
tion" and  dealing  with  the  extension  of  the  Hindu  method  of  evo- 
lution to  the  solution  of  numerical  equations  by  Vieta  and  Horner; 
and  a  more  extended  discussion  by  F.  Cajori  in  an  article  entitled 
"A  History  of  the  Arithmetical  Methods  of  Approximation  to  the 
Roots  of  Numerical  Equations  of  one  Unknown  Quantity,"  in  the 
Colorado  College  Publication  (Colorado  Springs,  Colorado)  for  1910. 
In  the  present  survey  we  purpose  to  study  the  origin  and  growth  of 
the  various  methods  by  bringing  out  the  illustrations  and  the  details 
of  explanation  as  far  as  possible  from  the  original  sources. 

The  writer  has  had  unusual  opportunities  for  studying  such  source 
materials  in  four  notable  libraries, — the  Columbia  University 
Library,  the  Teachers  College  Library,  the  New  York  City  Public 
Library,  and  the  private  library  of  Professor  David  Eugene  Smith, 
with  its  unique  collection  of  editions  of  the  Renaissance  and  early 
post-Renaissance  periods.  By  this  means  he  has  had  access  to  the 
early  editions  and  standard  translations  of  Euclid,  Archimedes, 
Heron,  Theon,  Brahmagupta,  Mahavira,  Bhaskara,  Omar  Khay- 
yam, as  well  as  the  latest  brochures  and  translations  of  Chinese  and 
Japanese  writers  on  mathematics.  He  has  also  been  able  to  consult 
such  works  as  those  of  Leonardo  of  Pisa,  Chuquet,  Pacioli,  Bom- 
belli,  Cardan,  Stevin,  Pitiscus,  Vieta,  Oughtred,  Hume,  Herigone, 
Wallis,  Newton,  Raphson,  Colson,  Rolle,  DeLagny,  Halley,  Taylor, 
Waring,  Euler,  Lagrange,  and  many  others  who  wrote  upon  the 
subject  previous  to  Horner. 


PRE-ALGORITHMIC  METHODS 

I.     EGYPTIAN  AND  BABYLONIAN  WORK  IN  FINDING  ROOTS 

We  formerly  thought  that  the  root  concept  originated  among 
the  Greeks,  but  recent  discoveries  show  us  that  the  rinding  of  roots 
of  numbers  and  the  solution  of  quadratic  equations,  and  even  the 
approximation  methods  employed  in  this  work,  go  as  far  back  as 
the  early  Egyptian  and  Babylonian  civilizations. 

On  the  Senkereh  tablets  (c.  2000  B.  c.)  are  tables  of  squares  and 
cubes,  which  shows  that  the  Babylonians  had  at  least  an  indirect 
notion  of  square  and  cube  roots.1  Three  quadratic  equations  are 
known  to  have  been  studied  in  this  period.  The  first  one  was  made 
known  when  Griffith  (1897)  published  the  mathematical  papyrus 
found  by  Petrie  in  Kahun.2  It  deals  with  areas  and  requires  the 
solution  of  the  equations  xy  =  12  and  x  :  y  =  I  :  }4,  stated,  of 
course,  in  rhetorical  form.  In  1900  Schach  discovered  in  a  Berlin 
papyrus  a  second  problem  in  quadratics,  requiring  the  solution  of 
the  equations  x2  +  y2  =  100  and  x  :  y  =  1  :  ^.  The  third  equation 
was  found  in  the  Kahun  papyrus  by  Schach  in  1903.  It  requires  the 
solution  of  the  equations  x2  +  y2  =  400  and  x  :  y  =  2  :  i}4.  The 
ancient  mathematician  solved  it  by  letting  x  =  2,  y  =  i}4\  this 
gives  x2  +  y2  =  b%\  since  \/6}4  =  2Y>  and  2x/2  =  }/%  of  20,  he 
found  that  x  =  2-8  =  16  and  y  =  i}4'%  =  12. 

The  same  method,  commonly  called  the  method  of  False  Position, 
is  also  used  in  solving  the  other  two  equations. 

No  irrational  roots  occur  in  these  equations.  But  there  are  indi- 
cations 3  that  the  Egyptians  had  a  definite  way  of  approaching  the 
square  root  of  non-square  numbers.  However,  they  seem  to  have 
been  unaware  of  their  irrational  quality. 

2.     EARLY     CONCEPTS     OF     HIGHER     ROOTS     AND     IRRATIONALS 
AMONG  THE  GREEKS 

That  some  roots  are  irrational  was  first  recognized  by  the  Greeks. 
From  measuring  the  areas  of  commensurable  squares  and  rectangles 

1  M.  Cantor,  Vorlesungen  iiber  Geschichte  der  Mathemalik,  Vol.  I  (hereafter  desig- 
nated Cantor,  I),  Leipzig,  1907,  p.  28. 

1  M.  Simon,  Geschichte  der  Mathemalik  im  Altertum,  Berlin,  1909,  pp.  41-42. 
8  Simon,  Geschichte,  pp.  43-53. 
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in  terms  of  unit  squares  they  proceeded  to  find  the  areas  of  non- 
rectangular  and  irregular  figures,  such  as  the  triangle  and  the  circle. 
This  they  did  by  finding  an  equivalent  square  (Quadrature).  Since 
they  neither  possessed  algebraic  symbols  nor  used  the  equational 
form  of  statement,  they  found  the  unknown  side  by  geometrical 
construction.  Some  of  the  sides  were  discovered  to  be  incommensu- 
rable quantities.  As  early  as  525  b.  c.  Pythagoras  recognized  the 
diagonal  length  of  the  unit  square  to  be  an  irrational  number. 
Theodorus  of  Cyrene  (c.  410  B.  C.)  and  Plato  {c.  380  b.  c.)  discussed 
the  irrational  nature  of  the  expressions  V3,  V5,  V7,  •  •  •  V17. 
In  Book  X  of  his  Elements  Euclid  (c.  300  b.  c.)  gave  a  comprehensive 
exposition  of  the  theory  of  irrationals.  Up  to  the  time  of  Archi- 
medes square  roots  were  almost  universally  evaluated  by  geometri- 
cal constructions;  by  him  they  were  often  evaluated  numerically 
by  a  process  of  approximation. 

The  Greeks  seem  also  to  have  been  the  first  people  to  conceive 
directly  the  notion  of  cube  root.4  This  also  arose  from  geometrical 
problems,  and,  like  the  square  root,  was  approached  geometrically. 
According  to  the  familiar  story  the  god's  demand  of  an  altar  twice 
as  large  as  the  old  one  necessitated  the  solution  of  the  equation 
x?  =  2as.  This  was  probably  solved  by  Menaechmus  (c.  365  b.  c.) 
by  the  aid  of  conic  sections.  Archimedes'  famous  problem  of  the 
section  of  the  sphere  led  to  an  equation  of  the  type  x3  —  ax'2  +  6=0, 
which  was  also  handled  geometrically.  Heron  of  Alexandria  (c. 
200  A.  D.)  was  the  first  one  to  make  a  practise  of  evaluating  cube 
roots  numerically;  this  he  did  in  the  course  of  his  practical  work  by 
a  process  of  approximation. 

3.    THE  APPROXIMATIONS  OF  ARCHIMEDES  AND  HERON  OF  ALEXANDRIA 

In  his  Measurement  of  the  Circle  Archimedes  (287-212  b.  c.) 
gives  among  others  the  following  approximations:5 

^  >  V3  >  —  ',    V349450  >  594;    V5472I32tV   >  2339-. 
780  153  8  4 

Nowhere  has  he  disclosed  his  process,  and  many  conjectures  and 
many  ingenious  restoration  formulas  have  been  advanced  as  ex- 

4  Tropfke,  Geschichle  der  Element ar-Mathemalik  I  (1),  Leipzig,  1902,  pp.  158-61; 
207-10;  223-24. 

'J  Cantor,  I  (1907),  p.  316. 
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planations.6  We  have  one  bit  of  information.  In  discussing  Heron 
of  Alexandria  (r.  200  a.  d.),  who  lived  four  centuries  later,  Eutocius 
tells  us  that  he  found  square  and  cube  roots  by  the  same  method  as 
Archimedes.  That  remark,  together  with  testing  for  certain  recon- 
struction formulas,  seems  to  indicate  that  at  times  he  used  methods 
of  Double  False  Position,7  as  in  evaluating  the  expression  y/2,  and 

at  other  times  he  used  the  formula  y/N  =  -  (a  +  — J,  as  in  find- 
ing the  square  root  of  3. 

Now,  we  have  Heron's  processes  described  in  his  own  words.  A 
tenth  century  manuscript  of  Heron's  Metrica  was  recently  discov- 
ered by  Dr.  R.  Schone  in  the  Serailo  Library  in  Constantinople.  In 
describing  how  to  find  the  square  root  of  720,  the  area  of  a  certain 
triangle,  Heron  says: 

"Now,  since  720  has  no  rational  square  root,  we  shall  find  the 
root  differing  from  it  by  a  very  small  error.  Since  the  next  larger 
square  to  720  is  729  and  its  side  is  27,  we  divide  720  by  27.  The  result 
is  26^.  To  this  add  27 ;  that  gives  53%.  One-half  of  this  is  26  Y  Yz . 
Therefore  the  next  root  of  720  is  26  Y*  Yi-  Multiplying  26  Y  Yi 
by  itself  gives  720V86,  so  that  the  error  amounts  to  only  Yse  of  unity. 
Should  we  wish,  however,  to  obtain  an  error  smaller  than  Mw,  we 
may  put  in  place  of  729,  the  value  we  have  just  found,  720  Vfee,  and 
if  we  do  this,  we  shall  find  an  error  much  smaller  than  Y>w." 

For  want  of  a  name  we  shall  call  this  the  "method  of  averages." 
If  a  is  a  rational  root  of  N,  then  N-t-a  =  a.  But  if  it  is  an  approxi- 
mation, and  a  little  too  large,  as  in  Heron's  example,  then 
N  -r-  a  =  au  a  little  smaller  than  the  exact  root.    Heron  takes  as 

the    final   approximation  ~(a  -f-  d\)  =      [a  +  —  J. 
2  2  V  a  / 

For  approximating  the  cube  root  *  Heron  used  the  method  of 

Double  False  Position.    He  explains  this  process  as  follows:  "We 

shall  now  explain  how  to  find  the  cube  root  of  100  units.    Take  the 

two  cubes  nearest  to  100,  one  larger,  the  other  smaller.    Note  how 

6  Cantor.  I  (1907),  pp.  315-17;  consult  S.  Giinther,  "Quadratische  Irrationaliteten 
der  Alten,"  in  Abhandl.  zur  Geschichle  der  Math.,  Vol.  IV. 

7  G.  Wertheim,  "Heron's  Ausziehung  der  irrationalen  Kubikwurzeln,"  in  Zeitschrift 
fur  Math.  u.  Phys.,  XLIX  (1899),  Hist.  Lit.  Abt..  pp.  1-3. 

8Curtze,  in  Zeitschrift  fiir  Math.  it.  Phys.,  XLII  (1897),  p.  117  119:  Wertheim,  in 
XLIV  (1899),  pp.  1-3;  Bibliotheca  Maihematica  (hereafter  called  Bibl.  Math.),  third 
series,  VIII  (1907-1908),  p.  412;  Cantor,  I  (1907).  P-  374- 
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much  larger  the  first  one  is,  namely  25,  and  how  much  smaller  the 
second  one  is,  namely  36.  Then  multiply  36  by  5;  the  result  is  180. 
To  this  add  100,  which  gives  280  [and  divide  180  by  280];  the  result 
is  %4.  Add  this  to  the  root  of  the  smaller  cube,  namely  4;  the  result 
is  4  %4.  Which  is  the  value  of  the  cube  root  of  100  units  as  nearly 
exact  as  possible." 

The  ambiguity  of  the  numbers  in  this  example  (5  maybe  -^125 
and  also  V125  -  100,  etc.)  excludes  a  precise  knowledge  of  Heron's 
process;  but  it  seems  to  differ  from  those  used  by  the  medieval 
writers.  Curtze,  Wertheim,  Cantor,  and  Enestrom  have  all  worked 
out  different  reconstruction  formulas  for  Heron's  process  for 
approximating  cube  roots. 

4.  theon's  method  of  exhaustion 
The  first  instance  in  extant  writing  of  a  method  of  finding  square 
root  similar  to  our  present  method  is  found  in  the  writings  of  Theon 
of  Alexandria  {c.  365  a.  d.).  He  gives  a  very  full  description  of 
how  to  evaluate  irrational  roots  in  sexagesimal  fractions.9  It  is  a 
geometrical  process  referred  to  a  figure.  In  it  he  applies  a  process 
of  exhaustion  to  Euclid's  geometrical  counterpart  of  the  expression 
(a  +  x)2  =  a2  +  2ax  +  x2.  This  is  followed  by  a  generalizing  rule. 
Expressed    in    modern    analytical    form,    if    a2  +  b  =  (a  +  x)2  = 

a2  +  2ax  -f-  x2  and  a  is  a  first  approximation,  then  x  =  —  approxi- 

2a 

mately,  and   Va2  +  b  =  a  -\ ,  approximately.     In  his  illustrative 

2a 

example  he  finds  the  square  root  of  45000  to  be  670  4'  45". 

Diophantus  did  brilliant  work  in  approximation,  but  this  was 

chiefly  in  connection  with  inequalities  and  indeterminate  equations, 

both  of  which  are  outside  of  our  topic. 

5.     SUMMARY 

1.  We  find  the  first  concept  of  roots  in  Egypt  or  in  Babylon;  in 
Egypt  we  also  find  the  solution  of  quadratic  equations  by  single 
false  position. 

2.  The  concepts  of  irrationality  and  of  higher  roots  are  first  found 
among  the  Greeks.  Theirs  are  also  the  earliest  records  of  a  scien- 
tific notion  of  limit  and  of  scientific  approximation. 

•  Gow,  History  of  Greek  Mathematics,  Cambridge,  1884,  p.  55. 
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ALGORITHMIC  METHODS  FOR  APPROXIMATING 

'I  HE  ROOT-,  OF   PURE  POWERS 

I.    THE  HINDU  METHOD  OF   EXHAUSTION 

The  superior  notation  and  convenient  symbolism  of  the  Hindus 
enabled  them  to  use  processes  of  inversion  to  a  degree  impossible 
to  I  neon.  In  faci  that  became  the  favorite  line  of  attack  in  much 
ol  their  mathematics.  Square  and  rube  roots  were  found  by 
bhaja ]  (b.  476  a.  d.)  from  the  inverse  of  the  formulas 
a1  4-  2<i\>  +  tr  and  <r  +  sa?b  +  3a/,2  +  h\ 

When  the  Hindu  decimal  notation  with  its  place  values  came 
into  use,  the  work  of  rinding  roots  was  so  arranged  that  an  applica- 
tion of  the  inversion  process  would  evolve  the  digits  of  the  root  in 
order  (whence  "evolution"  in  texts  on  arithmetic  and  algebra).  To 
eflfecl  this  the  number  was  divided  into  periods  and  the  work 
arranged  in  columns  and  lines.  This  "Hindu  method,"  as  it  was 
often  called  in  the  Middle  Ages,  is  described  for  us  by  Brah- 
magupta  (b.  598  a.  d.).2  Sridhara,3  in  his  Trisatika  (1040  a.  d.) 
gives  the  rules  in  full,  as  also  does  Bhaskara  (b.  11 14  a.  n.).4  The 
following  is  a  free  translation  from  Bhaskara,  by  Taylor: 

"The  firsl  place  on  the  right  is  called  ghana  or  cube;  the  two  next 
places  aghana  or  not-cube.  Subtract  the  cube  contained  in  the 
final  period  from  the  said  period;  put  down  the  root  of  the  cube 
in  a  separate  line,  and  after  multiplying  its  square  by  three,  divide 
the  antecedenl  figure  by  the  result,  and  write  down  the  quotient 
in  the  separate  line:  Then  multiply  the  square  of  the  quotient  by 
the  preceding  number  in  that  line  and  by  three,  and  after  sub- 
tra<  ting  the  producl  from  the  next  antecedent  figure  cube  the  said 
quotient,  and  subtract  the  result  from  the  next  antecedent  figure. 
Thus  repeat   1 1 n-  process  through  all  the  figures.    The  separate  line 

tains    t  lie    (    llbe    Root  ." 

1  <  antor,  I  ( [907  k  p.  625. 

II     11  oleb ke'a  translations  "I    Bialmiagupta  and    Bhaskara,  London.  1817. 

pp     ■,<)  81      Hereafter  referred  to  as  Colebrooke. 

1. .,•  e,  on  Tristika,  in  Bibl.  Math.  (3)  XIII  (1012   13).  P-  209. 
1      loi      translation  ol  Lilawali,  Bombay,  i8if>,  p.  20. 
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and  al-Qalasadi  (c.  1475). 10    This  rule  was  also  used  by  the  Jewish 
writer  Johannes  Hispalensis  (c.  1140).    Al-Karchi  (c.  1020)  used  a 

b 


different  formula,  namely  \/a2  -f-  b  =  a  +  •       How    he    de- 

2a  +  1 
rived  this  formula  is  not  known;  possibly  he  used  a  method  of  double 
false  position  and  interpolation  by  proportional  parts  similar  to 
the  one  used  later  by  Leonardo  of  Pisa  for  cube  roots.  Ibn  Albanna 
(c.  1300)  used  al-Karchi's  formula  for  b  >  a,  but  Theon's  formula 
for  b  ^  a.    Leonardo  used  Theon's  formula  with  a  corrective  sup- 


\2a/ 


plement,  as  follows:  \/a2  +  b  =  a  -\ —   .   A  formula 

2a         (      .    b  \ 

V  20/ 

frequently  used  in  the  Middle  Ages  was  y/A   =  -  VAa2.      In  the 

a 

Renaissance  period  Bombelli  (1572),  Cataldi  (1613),  and  Schwenter 

(1618)  approximated  irrational  square  roots  by  the  use  of  continued 

fractions. 

A  formula  for  approximating  irrational  cube  roots  was  invented 

by  Leonardo  of  Pisa.11   In  deriving  this  he  used  double  false  position 

and  interpolation  by  proportional  parts.    He  lets  a3  and  (a  -f-  i)3 

be  the  two  cubes  nearest  to  a3  +  b.    Then,  since  a3<    a3    -f-    b    < 

(a  +  i)3,  we  have  oy<  b  <  3a  (a  +  1)  +  1.    This  relation  gives  a 

criterion  for  the  proper  choice  of  a.    Also,  if  we  increase  a  by  unity 

we  increase  the  resulting  cube  by  3a  (a  -f-  1)  -f-  1.   What  increase  in 

a     will     increase     the     resulting     cube     by     b?      By    proportion 

he  obtains  as  an   answer    .      Therefore    \Va3  +  b  = 

3a(a+  1)  +  1 

b  12 

a  +  .       He  gives  these  illustrations  of  his  formula: 

3<z(a  +  1)  +  1 

10  Giinther,  Quadratische  Irrationatiteten  der  Allen,  pp.  45-46;  Cantor,  I  (iqo8), 
PP.  633,  641;  Hankel,  Zur  Geschichle  der  Mathematik  im  AUerthum  und  Mittelalter 
(hereafter  referred  to  aa  Hankel),  Leipzig,  1874,  p.  185. 

11  Boncompagni,  Leonardo  Pisano,  I,  pp.  370,  378,  380-81;  Cantor,  II  (1900),  pp. 
31-32;  Bibl.  Math.  (3)  II  (1902),  pp.  350-54- 

12  Juan  de  Ortega  (1512)   used  a  modification  of  this  formula,    viz.,    </ ai  4.  f,   = 

b 
a    +         ,      ,       '    See  Cantor,  II  (1900),  p.  388. 


v^oo  =  9  + 
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-1—  =  9  -  approximately,   and    -^2345  =   13  -f-  -*— 
271  3  547 


=  1 3-,  approximately. 
4 
To  find  the  roots  of  large  irrational  numbers  two  plans  were 
followed,  the  Hindu  arrangement  being  used  in  both  cases.  Leo- 
nardo illustrates  them  both.13    In  one  of  his  solutions,  V927435  = 


963  +  

321 


folk 


(")' 


2(962  +  II) 

\  321/ 


the  derivation  of  963  is  shown   as 


We  find  this  method  used  as  late  as  the  seventeenth  century 
in  the  Kholasat-al-Hisdb  14  of  Beha  Eddin  (c.  1600).  Another  plan 
is  followed  by  Leonardo  in  this  solution: 

/ l        /— l  o  l  o        l        l 

V7234  =  — v 72340000  = .  8505  -  =  85  -    — . 

100  1000  4  20  400 

This  plan  became  increasingly  popular  and  when  in  1539  a  man  of 
Cardan's  prestige  adopted  it  systematically  in  his  arithmetic,15 
accompanied  by  clear  and  concise  rules,  it  became  the  standard 
method  in  Europe. 

It  had  taken  a  thousand  years,  from  Aryabhata  to  Cardan,  to 
perfect  this  approximation  method.  But  it  is  the  most  effective 
and  most  widely  used  method  known.  The  methods  of  Vieta  and 
of  Horner  are  its  lineal  descendants. 

1S  Boncompagni,  Leonardo  Pisano,  I,  p.  355;  Cantor,  II  (1900).  p.  30;  Curtze. 
"Ueber  eine  Algorismus  Schrift  des  XII.  Jahrhunderts"  in  Abhandl.  Geschicht.  Math., 
VIII  (1898).  pp.  1-28. 

14  Taylor's  Lilawati,  Preface,  p.  14;  Hutton's  Tracts  II,  London.  1812,  p.  180.  A 
manuscript  copy  of  the  Kholasat-al-Hisdb  is  found  in  Professor  D.  E.  Smith's  private 
collection. 

15  Cardan's  Arithmetic,  chap.  23. 
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3.     SUMMARY 

1.  The  superior  symbolism  and  notation  of  the  Hindus  enabled 
them  to  invent  and  use  with  ease  formulas  based  on  the  principle 
of  inversion. 

2.  For  finding  the  rational  roots  of  numbers  too  large  for  hand- 
ling by  inspection  they  used  a  method  of  exhaustion.  It  was  based 
on  the  inversion  of  the  expressions  a2  -\-  2ab  -f  b2  =  (a  +  b)2  and 
a*  +  i,a2b  -\-  306-  +  b:i  =  (a  +  6)3,  and  on  the  decimal  notation  with 
its  place  values.  Its  essence  was  the  evolving  of  the  digits  of  the 
root  by  orders  (Evolution).  This  method  was  also  used  for  irra- 
tional roots  later  on,  and  European  scholars  standardized  an  adapta- 
tion of  it  for  approximating  irrational  roots  by  decimal  fractions. 

3.  The  five  principal  approximation  formulas  invented  and  used 
in  the  Middle  Ages  were: 

b  ,.,       ,^-r  b 


(a)    vV  +  6  =  fl+~;  (b)     vV  +  b  =  a  + 


2a  2a  +  1 


(c)    VA  =      y/Aa2\  (d)    vV  -f  b  =  a  + 


a  3a (a  +  1)  +  1 


(e)    </a*  +  b  =  a  + 


3fl(a  +  1) 


IV 
APPROXIMATIONS  FOR  SPECIFIC  PURPOSES 

I.    THE  RISE  OF  NUMERICAL  HIGHER  EQUATIONS  AMOM, 
THE  ARABS 

Although  the  Arabs  did  not  contribute  much  original  matter  to 
algebra  they  vitalized  it  and  enriched  its  contents  by  applying 
algebraic  operations  to  the  problems  of  Greek  geometry  and  to 
their  own  problems  in  astronomy  and  trigonometry.  This  led  them 
directly  to  numerical  higher  equations.  Archimedes's  problem  of 
the  section  of  the  sphere,  which  led  to  a  cubic,  was  first  attempted  by 
al-Mahani  (c.  860),  and  was  later  solved  successfully  by  al-Khazin 
(<~-  95°)  with  the  aid  of  conic  sections.  The  trisection  of  angles  and 
the  computation  of  sides  of  regular  polygons  led  to  other  cubics,' 
which  were  also  solved  by  conic  sections;  for  they  soon  came 
empirically  to  the  conclusion  that  cubics  could  not  be  solved  alge- 
braically. The  solution  of  numerical  cubic  equations  by  intersecting 
conies  was  the  greatest  original  contribution  to  algebra  made  by 
the  Arabs.  These  solutions  remained  unknown  to  the  Western 
world,  and  were  rediscovered  in  the  seventeenth  century  by  Des- 
cartes, Thomas  Baker  and  Edmund  Halley.  The  success  of  the 
Arab  scholars  in  this  field  may  have  deterred  them  from  trying 
methods  of  approximation.  What  they  might  have  done  in  this 
field  may  be  inferred  rather  than  judged  from  the  solitary  example 
left  us  in  Saracen  writings. 

2.    A  SOLUTION'  BY  AN  UNKNOWN  ARAB  SCHOLAR 

In  Miram  Chelebl's  annotations  of  Hugh  Beg's  astronomic! 
tables  (1498)  there  is  explained  a  method  of  solving  the  equations 
Px  =  x3  +  Q,  where  P  and  Q  are  positive  and  x3  is  small  in  com- 
parison with  Q}  Chelebi  explains  that  this  equation  was  used  for 
a  specific  purpose,  namely,  to  find  the  sine  of  i°.  With  the  Persian 
astronomers,  as  with  us,  the  ordinary  trigonometric  interpolation 

1  L'algebre  d'Omar  Alkayyami  (translated  and  edited  by  Woepcke),  Paris,  1851, 
PP-  54-57;  P-  82;  Hankel.  p.  277. 

'■  Cantor  I  (1907),  pp.  781-82;  Hankel,  pp.  289-93;  Braunmiihl,  Gesckickle  der  Vnr- 
lesungen  der  Trigonomrlrie  I,  Leipzig,  1900-3,  pp.  72-74. 


14  Approximating  Roots  of  Numerical  Higher  Equations 

would  fail  for  sin  i°.  The  solution  is  commonly  attributed  to 
Giyat  Eddin  al-Kaschi.  Braunmiihl,  however,  thinks  it  is  due  to 
al-Zarkali  (c.  1080). 

In  this  particular  equation  P  =  45'  and  Q  =  47r  Gp  8'  29",  where 
the  radius  length  ir  is  the  sexagesimal  unit  and  is  equal  to  60^.   The 
solution  is  as  follows: 
From  the  equation, 

x  = ;  for  a  relatively  small  x  we  have  x  =  -. 

P  P 

Then  x  =  a  +  r;  a  +  r  =  - =  a  -\ — 

P  P 

n      I     R   +   a" 

—  a  + . 


Suppose =  b,  with  remainder  S;  that  gives  R  =  bP  +  S 


T,  ...  .    ,    .  Q  +  (a  +  b  +  s)3 

Then  x  =  a  +  b  +  s;a  +  b  +  s= - 

P 

-  n    J_    R    +    (°    +    b    +    5)3    _  -1_    ^    +    ^    +    b    +    5)J 

—  u-  -j-  — —  a  -\-  — 


P 


P  P 

Suppose  the  last  fraction  equals  c ,  with  a  remainder  T. 
Then  x  =  a  +  £>  +  c-f--?,  and  so  on. 

The  actual  problem  is  to  evaluate  sin  i°  from  the  known  trigo- 
nometric relation  60^  sin  30  =  3.6or  sin  i°  —  4  sin3  i°.  Letting 
x  =  sin  i°  and  simplifying,  the  equation  takes  the  form  45^  x  = 
x3  +  y^  •  6or  sin  30.  The  Arabs  knew  the  last  term  to  equal 
4/  6»  8'  29"  =  Q. 

The  following3  is  Chelebi's  problem,  carried  out  to  seconds;  in 
the  Arab  text  it  is  carried  out  to  fourths. 

»  Hankel,  p.  292. 
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Divisor          P 

45r 

a 

b 

c 

Quotient        x 

I* 

2' 

49" 

Dividend       Q 

47r 

6P 

8' 

29" 

aP 

45 

R 

2 

A3 

i 

R     a" 

127 

bP 

90 

S 

37 

(a  +   &)3  -  a3 

6 

12 

(a  +  b)3  -  a3 

2234 

41 

cP 

2205 

T 

29 

That  the  method  was  effective  for  its  particular  purpose  is  shown 
by  Ulugh  Beg's  tables.  That  it  had  its  limitations  is  also  obvious: 
it  is  laborious;  and  it  is  restricted  to  equations  where  P  and  Q  are 
positive  and  x  is  relatively  small.  Vet  the  advance  in  the  theory 
of  equations  was  soon  to  remove  the  latter  limitations;  and  refine- 
ments and  modifications  by  men  like  Harriot,  Raphson,  and 
Halley  might  have  reduced  the  former.  Certainly  this  remarkable 
solution  was  the  feat  of  some  genius  seemingly  unsupported  and 
isolated. 

We  quote  Hankel's  appreciation: 

"This  beautiful  method  of  solving  numerical  equations  stands 
second  to  none  of  the  methods  of  approximation  invented  in  the 
Occident  since  the  time  of  Vieta  in  fineness  and  elegance.    Apart 
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from  the  solutions  of  square  and  cube  roots,  with  which  it  has  a 
similarity  of  principle,  it  is  the  first  method  of  successive  numerical 
approximations  which  we  meet  in  the  history  of  mathematics." 

3.    A  SOLUTION  BY  LEONARDO  OF  PISA 

The  last  statement  may  or  may  not  be  true.  There  were  at  least 
two  other  modes  of  approximation  in  the  thirteenth  century.  One 
was  a  well-developed  system  in  China  which  we  shall  describe  in 
the  chapter  on  Horner's  method.  The  other  was  the  method  em- 
ployed by  Leonardo  of  Pisa. 

At  a  scientific  tournament  in  Pisa  in  1225  he  was  asked  by  the 
Emperor  Frederick  II  to  solve  the  equation  jc3  -f-  2.x2  +  io.v  =  20 
by  Euclid's  method  of  ruler  and  compasses.4  After  several  unsuc- 
cessful attempts  he  changed  his  method  of  attack  and  showed  by 
rigorous  analysis  that  the  roots  of  this  equation  could  not  be  found 
by  ruler  and  compasses.  Then  he  finally  gave  this  approximate 
solution:  x  =  i°  22'  7"  42'"  33  1V  4  v  40  V1.  What  process  he  used  no 
one  knows.  But  after  seven  centuries  of  further  study  in  calcula- 
tion   we   differ    from    his    result  by  onlv    i^vl  or 

31,104,000,000 

As  with  the  solution  of  the  Arab  scholar,  so  with  this  one:  we  marvel 
alike  at  its  isolation  and  at  its  accuracy. 

Many  attempts  have  been  made  to  explain  Leonardo's  method. 
Zeuthen  thinks  he  used  a  method  of  false  position  with  interpolation 
similar  to  his  method  for  the  cube  root  of  numbers.  Others  think 
there  was  some  relation  between  his  work  and  that  of  the  Arab 
scholar.  Professor  David  Eugene  Smith  thinks  he  was  acquainted 
with  the  work  of  the  Chinese.  There  is  evidence  in  his  writings 
that  he  was  familiar  with  the  Chinese  method  of  approximating 
roots  of  indeterminate  equations.5 

4.    SUMMARY 

1.  The  Arabs  added  little  that  was  original  to  the  science  of 
algebra.  But  they  vitalized  it  by  applying  the  new  science  to 
problems  of  geometry,  astronomy,  and  trigonometry. 

*  Cantor,  II  (1900),  p.  46;  Boncompagni,  Intorno  ad  alcune  opere  di  L,  Fisano,  pp. 
17-19. 

5  Cantor,  II  (1900),  p.  26;  Curtze  in  Zeitschrift  filr  Math.  U.  Phys.  XLI,  Hist.  Lit. 
Abt.,  pp.  Si,  82. 
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2.  This  gave  risr  to  higher  numerical  equations,  many  of  which 
they  solved  by  intersecting  curves.  The  success  of  this  method 
undoubtedly  delayed  the  work  of  solving  them  algebraically  by 
approximation. 

3.  We  have  only  two  extant  algorithmic  solutions  by  approxima- 
tion from  this  period,  namely  that  described  by  Miram  Chelebi 
and  the  solution  of  Leonardo  of  Pisa.  They  were  both  effected  for 
a  particular  purpose. 


METHODS  INVOLVING  DOUBLE  FALSE  POSITION 
AND  PROCESSES  OF  INTERPOLATION 

I .  THE  GENERAL  NATURE  OF  THESE  METHODS 

The  first  successful  general  methods  of  approximating  roots  in 
complete  equations  were  those  involving  "double  false  position," 
and  the  history  of  these  methods  is  that  of  incomplete  equations 
repeated.  The  chief  problem  in  this  mode  of  attack  is  to  discover  a 
practical  way  of  finding  intermediate  ("mean")  values  that  will  con- 
verge to  the  true  root.  Three  chief  plans  sprang  up:  the  "averaging" 
rule  of  Chuquet;  the  methods  of  interpolation  by  proportional  parts 
used  by  Cardan,  Biirgi,  and  Pitiscus;  and  the  process  of  orderly 
evolution  of  digits  used  by  Stevin,  known  later  as  Stevin's  Rule. 

2.   chuquet's  rule  of  mean  numbers 

Chuquet's  Triparty  1  (1484)  was  little  known  outside  of  France. 
It  was  circulated  only  in  manuscript,  not  being  printed  until  1880. 
Yet  an  examination  of  the  textbooks  of  LaRoche,  Petro  Sanchez 
Ciruelo,2  and  other  writers  reveals  that  he  had  considerable  influence 
in  France.  His  was  a  mathematical  mind  of  a  high  order.  His 
symbolism  and  arrangement  of  work  were  much  ahead  of  his  time. 
In  the  theory  of  equations  he  held  advanced  views;  he  recognized 
negative  and  hinted  at  imaginary  roots,  and  gave  a  definite  state- 
ment for  the  law  of  signs. 

In  his  Rule  of  Mean  Numbers  he  states,  without  proof,  that 

ai  +  a2    .  a-\  en 

: — ; — r  lies  between  —  and  —     Also,  if  0i  and  02  have  the  same 

01    ~\-    02  01  02 

signs,  the  following  differences  will  have  the  same  sign: 

ax       d\  +  a*       a;02  +  a20i  ax  -f-  a2       #2       a\bi  +  a26i 

b~~    b,    +    02    =    0!(0!    +    02)    aild   61    +    02   "  h    =     02  (by    +    b2) 

The  method  is  perfectly  general,  and  he  applies  it  to  complete  and 

1  Triparty,  pp.  653-54  >n  Boncompagni,  XIII  (1880);  Cantor,  II  (1900),  pp.  351, 
359-60;  Boncompagni,  XIV  (1881),  pp.  416,  419. 

2  Cantor,  II  (1900),  p.  386. 
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incomplete  equations  alike.    An  integer  a  is  put  in  the  fractional 

a 
form-.    For  rational  roots  the  approximation  terminates  of  itself; 

for  irrational  roots  the  rapidly  increasing  denominator  insures  a 
more  rapid  convergence  than  the  arithmetic  mean  that  had  been 
used  by  others.3    As  an  illustration  we  may  take  his  solution  of  the 

equation  x2  +  x  =  39  — . 
hi 

In  this  case  we  see  that  x  =  -  is  too  small,  while  x  —  -  is  too  large. 
1  1 

Hence  the  first  mean  is  =  — .    By  substitution  this  proves 

1  -f-  1         2 

too  small.    Therefore  the  new  limits  are  —  and  -,  the  new  mean 

2  1 

being  — =  — .    This  second  mean  proves  too  small.    The  lim- 

2  +  1         3 

its  —  and  -  give  the  third  mean  — =  — .     This  is  also   too 

3  1  3  +  14 

23           6     .       2\  -f-  6        29 
small.     The  limits  —  and  -  give  — =  — ,  which   proves  too 

4  1  4+i  5 

2X  29  2\  4"  29  ^2 

large.    The  limits  —  and  —  give  the  mean  — =  — ,  and  this 

4  5  4+59 

proves  to  be  the  exact  root. 


3.  THE  REGULA  AUREA  OF  CARDAN 

The  first  general  method  of  approximating  roots  of  numerical 
higher  equations  to  be  printed  and  publicly  known  was  the  "Regula 
Aurea"  (Golden  Rule)  of  Cardan  (1501-1576)  published  in  his 
Ars  Magna  4  (1545),  six  years  after  the  publication  of  his  rules  for 
square  and  cube  roots.  It  is  built  on  the  basis  of  two  false  positions 
and  a  particular  mode  of  interpolation.  Though  he  uses  it  only  for 
equations  of  the  third  and  fourth  degree,  it  is  applicable  to  equa- 

'  An  adaptation,  by  no  means  an  improvement,  of  Chuquet's  method  was  worked 
out  by  DeLagny  in  his  Method  of  Mediation  (1697). 

4  Chapter  XXX,  entitled  "De  regula  aurea";  Cantor,  II  (1900),  pp.  505-507- 
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tions  of  every  degree.  He  restricts  the  rule  to  equations  where  the 
coefficients  are  positive  or  zero. 

The  following  is  his  derivation  of  the  rule,  presented  in  modern 
analytical  form: 

Let/(x)  =  k,  where /(x)  is  a  polynomial  in  x  arranged  in  descend- 
ing order.  Let  a  and  (c+  i)  be  two  positive  integers  such  that 
f(a)  =  k  —  b  and /(a  +  i)  =  k  +  b' .  Then  the  root  X  lies  between 
a  and  a+i.  That  is,  x  —  a  +  t(a  +  i  —  a)  or  x  =  (a  +  i)  — 
«(a  +  i  —  a),  where  the  corrections  t  and  e  are  positive  proper 
fractions.  The  first  may  be  called  the  additive,  the  second  the  sub- 
tractive  form.  Since /(a  +  i)  >  f(x)  >/(a),  we  have /(a  +  i)  — 
/'(a)  >  f(x)  —  /(a)  >  o   and   consequently 

f(x)-f(a)       _         k-(k-b)  b      <j 


f(a  +  i)  -  /(a)        (*  +&')_(*-&)       b  +  6' 

Take as  an  approximate  value  of  /.    Then  x  =  a  + 

&  +  &'  6  +  6' 

Since  <  I  our  next  additive  form  will  then  be  f  (  a  +  ) 

6  +  6'  V         b  +  b'J 

=  k  -  b". 

Placing  this  new  additive  value  in  the  subtractive  form  we  have 

,->  +  .)-.{  (.+  1)-(.  +  4-|p)}-.  +  .-6-^. 

In  the  subtractive  form  e  is  less  than  unity.    For  an  approximate 

.  /(a  +  !)-/(*)  (£  +  &')_£ 

value  set  e  = 


*.  +  ,)_/(.+    »  ) 


6  +  6'  ' 


-.    This  gives  the  third  approximation  x  =  a  +  I 


6'  +  b"  6  +  6' 

6' 


6'  +  V 


=  r,  say.    According  as  k  lies  between  /(/)  and  f{a)  or  be- 


tween f{r)  and  /(a  -f  i)  one  must,  in  continuing,  use  the  additive  or 
the  subtractive  process. 
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His   first   illustrative   example   is   the   solution   of   the   equation 
x4  +  3.r3  =  100,  the  root  of  which  he  finds  to  be  2 
lowing  is  his  solution  and  arrangement: 


£775 
4697' 


The    fol- 


Contrary  to  his  own  restrictions  on  the  signs  of  the  coefficients 
Cardan  uses  the  rule  in  the  equation  x3  =  6x  +  20,  with  disastrous 
results. 


4.    STEVIN  S  CONTRIBUTION' 

The  approximation  work  in  Stevin's  algebra  (1685)  might  with 
equal  propriety  be  discussed  in  connection  with  Vieta  or  Newton. 
For  Stevin  (1548-1620)  anticipated  Newton's  analytical  form  and 
treatment  as  well  as  Vieta's  idea  of  evolving  the  digits  of  the  root 
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of  a  complete  equation  order  by  order.  But  since  his  scheme  of 
deriving  the  digits  involved  the  successive  choice  of  two  positions, 
we  discuss  his  method  in  this  chapter.  We  can  best  explain  his 
process  by  recording  one  of  his  illustrations.5 

Suppose  x3  =  2>oox  +  33915024; 
if  x  =         i,  then         i3  <  300- 1         +  33915024; 
if  x  =       10,  then       io3  <  300- 10      +  33925024; 
if  x  =     100,  then     ioo3  <  300- 100    +  33915024; 
if  x  =  1000,  then  iooo3  >  300- 1000  +  33915024. 

Hence  x  lies  between  100  and  1000.  By  trying  x  =  200,  x  =  300, 
x  =  400,  he  finds  the  root  to  lie  between  300  and  400.  Similarly 
trying  310,  then  320,  then  330,  he  finds  x  to  lie  between  320  and 
330.  A  similar  procedure  for  the  units  gives  324  as  the  exact  value 
of  the  root.  Stevin  observes  that  an  irrational  root  can  be  approxi- 
mated to  within  any  desired  degree  of  accuracy  by  using  this 
method  and  his  scheme  of  computing  decimal  fractions.  This 
method  was  advocated  and  used  considerably  by  Albert  Girard  in 
his  Invention  Nouvelle  en  V Algebre  (1629).  Though  laborious,  the 
method  is  general.  "Stevin's  rule"  was  used  by  later  algebraists, 
like  Oughtred  (1631),  Kersey  (1673),  and  Saunderson,  in  connec- 
tion with  other  methods.6 

5.     THE  METHODS  OF  BURGI  AND  PITISCUS 

Pacioli's  Summa  (1494),  the  earliest  printed  work  on  arithmetic 
and  algebra,  gave  to  Renaissance  scholars  the  method  of  double 
false  position  and  interpolation  by  proportional  parts  used  by 
Leonardo  for  finding  roots  of  numbers.  Leonardo's  method  now 
came  to  be  employed  frequently  also  in  the  solution  of  higher 
numerical  equations,  and  Pacioli  even  used  it  in  an  exponential 
equation.7  This  method,  which  has  always  been  basic  in  finding 
trigonometric  functions,  was  used  by  Jost  Burgi  in  solving  equa- 
tions which  rose  out  of  his  work  in  trigonometry.8   We  give  his  solu- 

6  Stevin,  Les  CEuvres  Mathematiques  (editor  Girard,  1634),  I,  problem  77;  Cantor, 
II  (1900),  p.  628. 

e  N.  Saunderson,  Elements  of  Algebra,  II,  Cambridge,  1790,  p.  728. 

7  Pacioli,  Summa,  fol.  186  recto  to  188  recto;  Cantor,  II  (1900),  pp.  325-26. 

8  Cantor,  II  (1900),  pp.  626,  645-46,  648;  C.  I.  Gerhardt,  Geschichte  der  Math,  in 
Deutschland,  Munich,  1877,  pp.  76-82;  84-87. 
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tion  of  such  an  equation,  namely,  9 -303c2  -f-  27X4  —  gx9  +  x8  =  o. 

The  root  has  been  found  to  lie  between  0.68  and  0.69;  let  x  =  0.68 

+  h,  and  f(x)  =  o.    Then 

/(0.68)  =       0.0569                   /(0.68)  =          0.0569 
/(0.69)  =  -  0.0828  /(0.68  +  h)  =  o. 

Difference  0.0 1 : 0.1397  Difference  h:  o  .0569 

Hence  0.1397  :  0.0569  =  0.01  :h;    or  h  =  0.004,   and  x  =  0.684. 

Repeating    the    process    with    /(0.6840)    and   /(0.6841)    he    gets 

0.00140012  :  0.00056410  =  0.0001  :  h;  or  h  =  0.00004029. 

Hence  x  =  0.68414029. 

A  more  conventional  method  of  double  false  position  is  used  by 

Pitiscus    (1612),    a    pupil    of    Biirgi.      In    solving    the    equation9 

x3 
5,176,380  =  3.x  —  x3  he  sets  x  =  1,725,460  -\ .    Manifestly  x  > 

3 
1 ,725,460.   By  trial  he  finds  x  to  lie  between  1 ,730,000  and  1 ,740,000. 
Let  us  call  these  two  values  a\  and  a2  and  represent  the  equation  by 
fix)  =  k.   Then /(a,)  =  k  -  38,158  =  k  -  d,;/(a2)  =  k  -  9,061  - 
k  —  dz.   Then 

a-idi  —  a\ 


d\  —  a\ 


=  1.743. ii4- 


6.     SUMMARY 

1.  The  first  successful  general  method  of  solving  complete 
numerical  equations  by  approximation  was  given  by  Chuquet  in 
his  Triparty  (1484).  It  was  not  printed  until  1880  and  remained 
comparatively  little  known. 

2.  The  first  printed  method  was  Cardan's  Regula  Aurea  (1545). 
It  was  restricted  as  to  the  signs  of  the  coefficients  but  not  as  to  the 
degree  of  the  equation.    It  never  came  into  general  use. 

3.  In  1585  Stevin  explained  a  method  of  evolving  the  digits  of 
the  roots  by  order. 

4.  Pacioli  and  Biirgi  extended  Leonardo's  method  of  finding  cube 
roots  to  the  solution  of  complete  equations  of  higher  degree.  This 
method  and  the  conventional  method  of  double  false  position  ex- 
emplified by  Pitiscus  were  the  ones  most  generally  used  before  the 
year  1600. 

•  Pitiscus.  Trigonometriae,  pp.  50-53. 


VI 

VIETA'S  EXTENSION  OF  THE  HINDU  METHOD 
OF  APPROXIMATION  TO  COMPLETE  EQUATIONS 

I .   vieta's  method  of  approximation 

A  new  epoch  in  the  history  of  the  solution  of  numerical  equations 
begins  with  the  year  1600,  the  date  of  the  first  publication  cf 
Vieta's  De  Numerosa  Potestatum  Pur  arum  atgue  Adfectarum  ad  Exe- 
gesin  Resolutione  Tractatus.  In  it  was  explained  a  systematic  and 
concise  way  of  solving  any  equation,  complete  or  incomplete,  by 
methodically  finding  the  successive  digits  of  the  root,  beginning 
with  the  highest  order.  Why  no  one  before  Vieta  should  have 
thought  of  applying  to  the  solution  of  complete  equations  the 
Hindu  method  of  finding  roots  of  large  numbers  may  seem  strange 
to  us  unless  we  reflect  that  useful  inventions  usually  seem  very 
natural  after  they  are  invented.  At  any  rate,  this  method  had  been 
applied  to  large  numbers  for  a  thousand  years,  had  appeared  in 
definite  forms  for  four  hundred  years,  and  had  been  perfected  in  its 
finest  details,  before  any  one  thought  of  extending  it  to  complete 
equations.  It  was  the  first  comprehensive  method  of  solving  such 
equations  that  had  been  attempted,  and  it  involved  no  restrictions 
as  to  terms,  signs,  or  degree,  for  it  applied  to  a  form  into  which 
any  equation  can  be  changed.  Vieta  calls  an  equation  "duly  pre- 
pared" if  it  has  the  form  x"  -f-  ax  _I  +  bx"~2  -f   =   K  where 

K  is  positive  and  all  the  coefficients  are   integers. 

De  Numerosa  Potestatum  was  first  printed  in  Paris  in  1600 
for  private  circulation.  In  1646  Franz  van  Schooten,  of  Leyden, 
again  printed  it  in  a  one-volume  collection  x  of  Vieta's  works. 
It  contains  sixty-six  pages,  of  which  the  first  ten  deal  with  in- 
complete equations.  He  arranges  the  solution  in  a  way  slightly 
different  from  the  conventional  one  in  order  the  better  to  show 
the   analogy   between    the   pure    and    the    affected    powers.2     We 

1  Vielae  Opera  Mathematica,  Leyden,  1646;  hereafter  referred  to  as  Vieta. 

2  Vieta  classified  powers  and  equations  as  pure  and  affected  ("ad-fected").  The 
names  are  geometrical  in  origin;  for  Vieta  considered  all  the  terms  in  an  equation  homo- 
dimensional.  Thus  in  the  "affected"  equation  x3  +  30X  =  K,  the  "resolvend"  K  is 
considered  three-dimensional;  and  30  is  the  "plane  coefficient"  in  the  "adjoined  solid" 
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insert  the  solution  for  one  pure  equation  and  two  types  of  affected 
equations,"  viz. : 

(1)  xz  =  157464;  where  x  =  54. 

(2)  .v3  +  30*  =  14,356,197;  where  x  =  243. 

(3)  95,400.1;  +  x3  =  1,819,459;  where  x  =  19. 

The  arrangement  is  that  of  Vieta.  The  explanations,  too,  are  a 
close  paraphrase  in  modern  notation  of  his  marginal  and  other 
explanations  rhetorically  expressed. 

(1)    Let/(*)  =  Kiorx*  =  157,464; 

if  x  =  a  +  b,  then  a3  +  i,a2b  4-  $ab2  +  b*  =  157,464. 
"Proponantur  1  C,  aequari  157,464. 

I.    Eductio  lateris  singularis  primi 


K  :  Resolvend 

a3 

K  -f(a):  Residual 


157 

464 

125 

32 

464 

(a 


32 

464 

7 

5 

15 

7 

65 

30 

2 

40 

64 

32 

464 

II.    Eductio  lateris  singularis  secundi 
K-f(a):  Resolvend 
3a2 
3o 

Divisor  7  65 

3a  "b 
3ab2 
b* 
Subtrahend 

Itaque  1  C,  aequetur  157,464  fit  1  N  54." 

(2)    Let/(x)  =  K  for  x3  4-  30*  =  14,356,107; 

if  *  =  a  4-  b,  then  (a3  4-  $a2b  +  3«62  4-  b*)  +  3o(a  4-  b)   = 
14,356,197. 

"Proponantur  1  C  -\-  30  N,  aequari  14,356,197. 

30A.-.  The  names  "pure"  and  "affected"  continued  in  use  until  the  latter  part  of  the  19th 
century.  Since,  in  the  modern  theory  of  equations,  the  basic  form  is  aox"  +  aix"  ~  l  + 
...  +  K  =  o,  recent  writers  use  complete  and  incomplete  as  being  more  descriptive. 
In  discussing  the  works  of  Vieta  and  his  immediate  successors,  we  shall  employ  the 
terminology  used  by  them. 

3  Vieta.  pp.  167-68;  177-78;  179- 
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I.    Eductio  lateris  singularis  primi 
Plane  coefficient 
A":  Resolvend  14  356  197  (a 

a* 

30a 

f(a) 

K-f(a):  Residual 


3 

0 

14 

356 

197 

8 

6 

0 

8 

006 

0 

1   6 

350 

197 

I.    Eductio  lateris  singularis  secundi 


Plane  coefficient 

K  -f(a)  :  Resolvend     6 

3«2 

3a 

Divisor 

&b 

2,ab' 

b> 

30b 

Subtrahend 

K-f(ab) 


3°. 

\     6 

35o 

197 

1 

2 
6 

1 

260 

30 

4 

8 

96 
64 
1 

20 

5 

825 

20 

524 

997 

<b 


III.    Eductio  lateris  singularis  tertii 


Plane  coefficient 
Resolvend 

3^2 
3<* 

Divisor 
$a*b 

tab* 

b* 

306 

Subtrahend 


3? 

524 

997 

172 

8 

72 

173 

550 

518 

4 

6 

48 

27 

90 

524 

997 

(b 


Itaque  1  C30  N,  aequetur  14,356,197  fit  1  N  243. 
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(3)    Let/(x)  =  K  for  95,400*  +  x*  =  1,819,459; 

if  x  =  a  +  b,  then  95  40o(a  +  b)  +  (a3  +  2>a}b  +  3a*2  +  &3) 
=  1,819,459. 
"Paradigma  cum  solidum  adfectionis  sub  latere  majus 

est  cubo. 

I.    Eductio  lateris  primi  inanis  ante  devolutionem 
Plane  coefficient 
K:  Resolvend 

Quoniam  9  major  est  unitate,  fit  devolutio. 


9 

540 

0 

1 

819 

459 

II 


Eductio  lateris  singularis  primi  post  devolutionem 


954 

00 

I 

819 

459 

954 

1 

00 

955 

864 

459 

Divisor 

K:  Resolvend 

95  400a 
a8 

/(a) 

K-f(a):  Residual 

III.    Eductio  lateris  singularis  secundi 
Plane  coefficient 
K~f(a):  Resolvend 

3«2 

Divisor 
95  400b 
2>a*b 
3ab> 
b* 
Subtrahend 

Itaque  si  95,400  N  +  iC,  aequentur  I 
1  N  19." 


(a 


95 

400 

864 

489 

3 

3 

95 

730 

858 

600 

2 

7 

2 

43 

729 

864 

459 

(*-9 


119,459     fit 


In  a  similar  fashion  he  illustrates  the  solution  of  equations  of  the 
fourth,  fifth,  and  sixth  degrees.  Vieta's  work  involves  mainly  four 
points : 
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1.  The  first  step  is  to  ascertain  the  number  of  digits  in  the  root. 
For  pure  powers,  of  course,  this  is  easy.  For  affected  powers  there 
will  also  be  as  many  digits  in  the  root  as  there  are  periods  in  the 
resolvend  if  the  highest  power  is  dominant;  that  is,  if  the  pure 
power  is  larger  than  the  adjoined  solid  or  solids.  This  will  generally 
be  indicated  by  the  size  of  the  coefficients.  Thus  examples  (i)  and 
(2)  will  have  roots  of  3  digits.  A  more  difficult  situation  arises  if 
the  highest  power  is  not  dominant.  Thus  in  example  (3)  the  re- 
solvend has  3  periods,  yet  the  root  has  only  2  digits.  Such  a  case 
he  calls  "devolution."  If  the  cube  and  the  adjoined  solid  had  been 
of  opposite  signs  there  would  have  been  more  digits  in  the  root  than 
periods  in  the  resolvend.  Such  cases  are  called  "anticipations." 
Vieta  and  all  his  followers  found  these  cases  full  of  practical  diffi- 
culties and  he  offers  many  suggestions  for  reducing  and  circum- 
venting them.    One  device  is  shown  in  our  third  illustration. 

2.  With  Vieta  the  fundamental  operation  for  finding  the  roots 
of  pure  as  well  as  affected  equations  is  the  expansion  of  the  powers 
of  a  binomial.4 

3.  For  each  approximation  a,  Vieta  takes  the  residual  K-f(a)  as 
the  new  resolvend.    The  entire  residual  is  put  down  each  time. 

4.  Vieta  gets  his  trial  divisors  by  taking  the  sum  of  all  the  co- 
efficients of  the  powers  of  b  except  the  coefficient  of  b",  which  is 
invariably  unity.  Both  Cantor  and  Hankel  have  given  erroneous 
explanations  of  Vieta's  divisors.  The  mistake  seems  to  come  from 
representing  Vieta's  process  in  the  compact  equational  form  of 
Newton's  scheme.  At  least  they  attribute  to  him  Newton's  divisor. 
Following  Hankel's  and  Cantor's  explanation,  the  first  divisor  in 
Example  (1)  should  be  120,030.  Vieta's  divisor,  however,  is 
120, 630.5  In  this  connection  it  should  be  said  that  both  these  his- 
torians give  an  erroneous  impression  of  the  plan  of  Vieta's  method. 
They  lead  one  to  suppose  that  Vieta's  method  of  evolving  digits 
is  the  same  as  Newton's  method  of  substituting  values  in  equa- 
tions of  series;  and  Hankel  says  as  much.  But  the  two  differ  in 
origin  and  form,  and  in  certain  of  the  steps. 

In   developing   his  approximation   method   Vieta   uses   for  con- 

4  Cantor's  exposition  of  Vieta's  process,  II  (1900),  pp.  640-41,  where  he  sets  x  = 
x\  +  Xt  +  x3.  is  apt  to  be  misleading.  Vieta  did  not  use  that  substitution,  though  his 
expositors  Hume  (1636)  and  Herigone  (1642)  used  it  occasionally. 

6  Hankel,  p.  370;  Cantor,  II  (1900),  p.  641. 
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venience  equations  having  rational  roots.  But  he  adds6  that  irra- 
tional roots  can  be  approximated  in  the  same  way  after  first  trans- 
forming the  equation  into  one  whose  roots  are  decimal  multiples 
of  the  roots  of  the  given  equation.  This,  of  course,  is  nothing  but 
an  extension  of  Cardan's  rule  for  irrational  numbers.  Thus,  in  the 
equation  x*  =  2,  the  root  is  found  from  solving  x3  =  2,000,000, 
;  likewise  if  .r3  +  6oo.r  =  8,000,  the  root  correct  to  hun- 
dredths is  found  by  solving  .vn  -+-  6o,ooo.r  =  8,000,000. 

2.     THE  MODIFICATIONS  OF  VTETA's  METHOD  BY  HARRIOT 

It  was  not  in  his  homeland  that  Vieta's  method  met  with  the 
greatest  appreciation.  The  people  of  the  British  Isles  have  always 
been  partial  to  mathematical  processes  dealing  with  numerical 
calculation.  Three  English  scholars,  Thomas  Harriot  (1560-1621), 
William  Oughtred  (1574-1660),  and  John  Wallis  (1616-1703),  all 
noted  as  teachers  and  writers,  became  the  chief  protagonists  of  the 
new  method. 

In  his  Artis  Analyticae  Praxis  7  (written  before  1621  and  pub- 
lished posthumously  in  1631)  Thomas  Harriot  introduces  two  devices 
for  expediting  the  calculation  in  Vieta's  method.  One  of  these  is 
the  elimination  of  the  second  term;  the  second  is  his  Canones 
Directorii,  a  collection  of  the  different  cases  of  numerical  solutions, 
ready  and  arranged  for  use  under  various  canonical  equations, 
showing  the  elements  necessary  to  form  the  several  subtrahends  and 
resol vends.  This  plan  gained  much  vogue  and  we  find  Raphson, 
Halley,  and  DeLagny  doing  the  same  thing  in  connection  with 
Newton's  method. 

In  the  actual  computation  he  uses  Vieta's  method  unchanged. 
We  find  no  basis  for  De  Morgan's  statement  that  he  formed  "only 
so  much  of  the  divisor  as  is  necessary  for  the  determination  of  the 
next  figure."  Nor  does  Wallis,  who  is  very  partial  to  Harriot,  men- 
tion that  fact.  We  subjoin  the  first  two  steps  of  an  example  from 
page  156  of  the  Praxis,  which  is  quoted  in  full  by  De  Morgan. 

6  Vieta,  p.  228. 

7  The  writer  has  been  unable  to  consult  the  original  and  depends  for  his  description 
on  the  following  secondary  sources:  De  Morgan's  article,  "Notices  on  the  Progress  of 
the  Problem  of  Evolution,"  in  Companion  to  the  British  Almanac.  1830,  pp.  34  52i 
Hutton's  Tracts,  II,  pp.  278-86;  Cantor,  II  (1900);  Henry  Stevens.  Thomas  Hariot  and 
his  Associates,  London.  [900;  J.  Wallis,  A  Treatise  of  Algebra  both  Historical  and  Prac- 
tical, London, 1685. 
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"Sit  aequatio  numerose  proposita 

aaaa  -  10240a  +  6254a  =  19633735875 

Radix  universalis  successive  educenda. 


Homogeneum  resolvendum 


hhhh 


3    7   5 
19633735875 


ggg 

6254 

-If 

—  1024 

bbb 

-  27 

+ 

27006254 

Divisor  A 

26903854 

gggb 

18762 

Radix  sing,  prima  b  =  3, 

-ffbb 

—  9216 

bbbb 

81 

+ 

81018762 

Ablatitium     Ab 

80097162 

Radix  singularis 

3 

Homogeneum  residuum  resolvendum 

1 1 62401 9675 

ggg 

6254 

-IJ 

—  1024 

-2ffb 

—  61440 

Radix  sing,  decupl.  b  =  30, 

4bbb 

108000 

6bb 

5400 

46 

120 

+ 

1 13526254 

- 

62464 

Divisor  B 

112901614 

gggc 

43778 

-ffcc 

-  50176 

-2ffb€ 

—  430080 

Radix  sing,  secund.  c  =  7 

^bbcc 

756000 

6bbcc 

264600 

^bccc 

41 160 

cccc 

+ 

2401 

1064204778 

- 

480256 

Ablatitium    Be       10594022 1 
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Harriot's  canonical  form  is  aaaa  —  ffaa  +  ggga  =  hhhh.  He  lets 
a  =  b  +  c.  Then  the  transformed  equation  becomes  (bbbb  + 
4bbbc  -f  6bbcc  +  \bccc  -f  cccc)  —  //  (bb  +  26c  +  cc)  +  ggg  (6  +  c) 
=  hhhh.  By  Vieta's  method  the  divisor  would  be  (ggg  —  ff  — 
2ffb  +  ^bbb  +  6bb  +  46) ;  and  this  is  identically  the  divisor  in  the 
illustration  quoted,  both  in  the  second  and  the  third  step. 

Harriot's  Praxis  was  the  most  analytical  algebra  that  had  been 
written  up  to  that  time.  His  lucid  treatment  of  the  theory  of  equa- 
tions, his  improved  symbolism,  and  the  inspiration  of  his  person- 
ality to  his  students  were  accountable  in  no  small  measure  for  the 
popularity  of  Vieta's  method  in  England,  and  it  was  often  called 
after  him  "Harriot's  method." 


3.     THE  IMPROVEMENTS  OF  OUGHTRED  AND  WALLIS 

No  one  did  more  to  popularize  the  new  method  than  did  the 
clergyman  mathematician  William  Oughtred.  This  he  accomplished 
by  giving  private  tuition  to  ambitious  young  men  and  these  spread 
his  teachings  throughout  Great  Britain;  among  them  were  Seth 
Ward,  Christopher  Wren,  and  John  Wallis.  An  even  more  potent 
influence  was  his  textbook,  the  Clavis  Mathematica.  Before  the 
year  1700  five  editions  of  this  little  volume  had  gone  through  the 
press. 

The  first  edition  appeared  in  1631,  the  same  year  in  which 
Harriot's  Praxis  was  published.  The  Clavis  appeared  first.  In  it 
there  is  no  reference  to  Harriot,  nor  any  use  made  of  his  improve- 
ments: the  two  men  evidently  knew  nothing  of  one  another's  work. 
Oughtred  uses  Vieta's  method  throughout.  In  "locating"  the  root 
he  depresses  the  equation  and  uses  a  method  similar  to  that  of 
Stevin  for  obtaining  the  first  digits  of  the  root.  It  is  interesting 
to  observe  that  he  uses  the  decimal  bar  and  logarithms  in  his 
solutions. 

John  Wallis  made  this  method  more  practical  by  arranging  the 
work  in  a  compact,  yet  clear  form,  so  as  to  show  the  steps  in  the 
computation.  He  differed  from  Vieta,  Harriot,  and  Oughtred  in 
not  rejecting  the  coefficient  unity  in  his  divisor.  To  help  discover 
cases  of  devolution  and  anticipation  he  employs  two  sets  of  periods, 
one  for  the  resolvend  and  another  for  the  coefficient;  at  times  he 
uses  Stevin 's  rule  for  this  purpose. 
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4.     THE  CONTRIBUTION  OF  HUME 

In  France  the  method  of  Vieta  made  slower  progress.  In  1636, 
James  Hume,  an  emigrated  Scotchman,  published  in  Paris  a  work 
with  the  title  Algebre  de  Viete,  d'une  methode  nouvelle,  claire,  et  facile. 
It  is  an  abridged  presentation  of  Vieta's  analytical  works,  repro- 
ducing both  his  methods  and  his  illustrations.  However,  he  intro- 
duces one  modification :  he  uses  only  so  much  of  the  divisor  and  of 
the  dividend  as  is  needed  to  compute  the  next  digit.  Thus  in  the 
example  (2)  which  we  quoted  from  Vieta  he  uses  12  for  the  divisor 
and  6350  for  the  dividend,  where  Vieta  uses  126,030  and  6,350,1 97. s 

Vieta's  method  was  adopted  by  P.  Herigone9  in  his  textbook 
Cursus  Mathematicus,  published  about  this  time;  but  neither  Hume's 
nor  Harriot's  improvements  were  used. 

5.    SUMMARY 

1.  The  Hindu  method  of  finding  the  roots  of  large  numbers  was 
extended  by  Vieta  to  the  solution  of  complete  equations  in  his  De 
Numerosa  Potestatum  (1600). 

2.  The  method  was  first  taken  up  in  England  by  Harriot  (some 
time  before  1621),  who  modified  it  by  using  canonical  equations  in 
which  the  second  term  had  been  eliminated. 

3.  Oughtred  (1631)  and  Wallis  (1685)  improved  the  arrange- 
ment, but  otherwise  used  Vieta's  method  unmodified,  except  that 
Wallis  included  the  coefficient  unity  in  his  divisor. 

4.  In  France  James  Hume  (1636)  used  only  so  much  of  the 
dividend  and  divisor  as  was  necessary  to  evolve  the  next  digit  of 
the  root.  Herigone's  algebra  (1642)  was  the  first  French  text- 
book to  popularize  Vieta's  method  in  France. 

8  Hume,  Algebre  de  Viete,  pp.  234-38. 

9  Herigone,  Cursus  Mathematicus,  II,  Paris,  1642,  pp.  270-76. 


VII 
NEWTON'S  METHOD  OF  APPROXIMATION 

I .  ADVANCE  IN  THE  THEORY  OF  EQUATIONS 

The  seventeenth  century  saw  a  remarkable  advance  in  the 
theory  of  equations.  The  symbolic  notation  introduced  by  Vieta  and 
improved  by  succeeding  algebraists,  notably  by  Harriot  and  Des- 
cartes, made  the  symbolic  equation  a  chief  instrument  of  analysis. 

When  Harriot  put  the  equation  in  the  form  a0xn~l  +  a,\Xn~l  +   

+  K  =  o,  the  eye  began  to  aid  the  mind  in  grasping  the  relations 
of  the  parts  of  an  equation.  The  law  of  signs,  used  by  Harriot, 
enunciated  by  Descartes,  and  refined  by  Newton,  gave  an  effective, 
though  limited,  criterion  for  the  number  and  nature  of  the  roots 
in  an  equation.  The  invention  of  coordinate  geometry  by  Descartes 
in  1637  and  the  independent  discoveries  of  the  differential  calculus 
by  Newton  and  Leibniz  a  few  years  later  gave  a  powerful  impetus 
to  work  in  theory  of  equations  and  in  practical  numerical  calcu- 
lation. 

For  purposes  of  numerical  approximation  one  of  the  tools  made 
most  effective  by  the  improved  symbolism  was  the  infinite  series 
in  an  equation.  Series,  as  a  means  of  exhaustion  and  of  summation, 
had  been  successfully  used  for  approximating  numerical  values 
since  the  time  of  Archimedes.  But  with  the  introduction  of  sym- 
bolic algebra,  series  gained  an  effectiveness  undreamed  of  by  the 
ancients.  Pioneers  in  the  use  of  equations  in  series  for  purposes  of 
numerical  approximation  were  Wallis  (1656),  Lord  Brouncker 
(1620-1684),  Gregory  (1667),  and  Newton  (1669). 

2.     A  DESCRIPTION  OF  NEWTON'S  PROCESS 

It  was  while  finding  the  area  under  curves  by  his  newly-invented 
"method  of  fluxions"  that  Isaac  Newton  (1 642-1 727)  discovered  a 
new  and  most  necessary  application  of  series.  Following  the  lead 
of  Wallis  he  transformed  such  fractional  and  radical  expression  as 


y  =  and  y  =  \/a?  +  x2  into  integrable  form  by  expanding 

b  +  x 
them  into  equations  in  series  by  ordinary  division  and  finding  of 
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roots.  This  scheme  he  extended  to  affected  equations.  "Since  the 
whole  difficulty  lies  in  the  solution,"  he  says,  in  explaining  it,  "I 
shall  first  illustrate  the  method  I  use,  in  a  numerical  equation."  l 

Such  then  is  the  historical  setting  of  the  beginnings  of  a  type  of 
approximation  by  successive  substitutions  that  is  almost  universally 
known  as  "Newton's  method  of  approximation."  It  was  entirely 
incidental  to  his  work  in  quadrature  by  calculus.  The  originator  of 
this  famous  method  gives  only  a  single  numerical  example,  and  it 
is  illustrative  and  introductory  to  the  more  general  problem  of 
deriving  a  method  for  finding  the  root  of  literal  equations  by  means 
of  equations  in  series.  Some  historians  notwithstanding,  it  is 
hardly  correct  to  say  that  it  was  a  further  development  or  an  adap- 
tation of  Vieta's  method.  Of  course,  Newton's  method  includes 
Vieta's  ideas  of  limit,  substitution,  and  exhaustion;  but  these  seem 
inherent  to  any  method  not  using  double  false  position.  In  some 
of  its  features  Newton's  solution  resembles  the  approximation 
processes  in  Chelebi's  annotations  and  of  Stevin.  Later  algebraists2 
sometimes  used  "Stevin's  rule"  of  substituting  1,10,100  etc.,  to  find 
the  integral  part  of  the  root.  Newton  himself  says  of  this  method: 
"Whether  this  method  of  solving  equations  be  commonly  practised 
or  not,  I  cannot  tell;  but  surely  to  me  it  appears  simple  in  com- 
parison with  other  methods  and  more  suited  to  practise."3 

Newton's  solitary  numerical  example,  the  solution  of  yz  —  iy  —  5 
=  o,  has  become  classic.  Lagrange,  RjUffini,  and  Horner,  all  demon- 
strated the  virtue  of  their  methods  by  applying  it  to  Newton's 
equation,  and   it  has  become  a   favorite  problem   in   text-books. 

This  process  of  approximating  by  successive  substitution  in  equa- 
tions is  thus  explained  by  Newton: 

Suppose  that  an  approximate  value  y  =  2  has  already  been 
found,  differing  by  not  more  than  one  tenth  of  itself  from  the  true 
root.  Then  y  =  2  +  p.  Substituting  this  value  in  the  original 
equation  gives  p3  +  6p2  +  lop  —  1  =  o.  Since  p  is  a  fraction  we 
can  reject  all  but  the  linear  terms,  or  lop  —1=0  approximately. 
That  is  p  =  0.1,  approximately;  or  p  =  0.1  +  q.  Continuing  this 
process  as  far  as  desired  we  get  several  values  q,  r,  s,  etc.    And 

1  Newton,  Analysis  by  Equations  of  an  infinite  Number  of  Terms,  translated  and 
edited  with  commentaries,  by  John  Stewart,  London.  1745,  pp.  321-28,  357. 
i  N.  Saunderson,  Elements  of  Algebra,  II,  Cambridge,  1790,  p.  728. 
3  Newton,  Analysis  by  Equations,  p.  330. 
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y  —  £  +  5  +  r  +  5+   If  at  any  time  the  coefficients  are 

such  that  one  doubts  the  accuracy  given  by  the  linear  terms,  one 
may  advantageously  use  the  quadratic  part  of  the  equation.  Thus 
6p'2  +  lop  —  I  =  o  would  give  a  more  accurate  value  of  p  than 
lop  —  i  =  o.  In  such  cases,  use  the  lesser  root  of  the  quadratic. 
Newton  sets  up  his  solution  in  paradigm  form,  as  follows:4 


•y3  — 

-2)>-5=o 

+  2,10000000 
-0,00544853 
+  2,09455147  =? 

2+p  =  y 

+y3 

+  2y 

-5 

sum 

+  8+i2£+6£2+/>3 
-4    -2P 

-5 

-i  +  iop+6p2+p3 

o,  i+q  =  p 

+ps 
+  6£2 
+  iop 

—  1 

+0,001  +o,03g+o,352+g3 
+0,06  +1,2     +6,0 
+  1,       +10 

sum 

+0,061  + 1 1 ,23g  +  6,3g2+g3 

-0,0054  +  ^  =  9 

+  6,392 

+  11,239 

+0,061 

sum 

+0,0001 83708  —  o,o68o4r+6,3r2 
—  0,060642       +11,23 
+0,061 

+0,000541708  +  1 1 ,16196^+6,3^ 

—  0,00004854  +  5  =  r 

He  next  takes  up  the  more  general  method  of  literal  equations.8 
"These  things  being  thus  shown  in  numbers,  let  y3  +  a2y  —  2a3  + 
axy  —  x3  =  o  be  proposed  to  be  solved."    Using  the  same  method, 

setting  v  =  a  +  p, he  derives  p= +  q;  q  =  ^    +  r,  etc.,  and 

4  64a 

y  =  a_*+il_  +  111*3  +   509*4  + A  practical  difficulty 


4  Newton,  Analysis  by  Equations,  pp.  328-329. 
s  Newton,  Analysis  by  Equations,  pp.  330-.343- 
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in  getting  the  supplementary  approximation  is  encountered  when 
the  coefficients  are  letters.  He  reduces  this  by  arranging  the  terms 
according  to  their  dimensionality;  this  arrangement  has  been  called 
Newton's  Parallelogram.  The  value  of  y  will  converge  when  coeffi- 
cients are  so  chosen  that  x  is  small;  and  the  smaller  they  are,  the 
more  rapidly  the  series  converges.  Later  algebraists  sometimes 
found  numerical  approximations  by  using  as  forms  such  solutions  of 
literal  equations.  But  Newton,  as  might  be  expected,  did  not  sub- 
stitute numerical  values  of  x  until  after  the  integration.  With  him 
the  objective  was  the  area  of  the  curve  and  not  the  numerical 
value  of  y. 

Newton's  description  of  his  method  was  contained  in  an  unpub- 
lished tract  entitled  De  Analyst  per  aequationis  numero  terminorum 
infinitas  which  was  shown  to  Barrow  and  Collins  in  1669.6  His 
method  was  first  given  to  the  public  in  the  1685  and  1693  editions 
of  Wallis's  algebra.  In  1704  and  171 1  it  was  published  in  Newton's 
own  works  in  Latin  and  in  1736  and  1745  it  appeared  in  Colson's 
and  Stewart's  translations.  Newton's  name  was  always  linked  with 
his  particular  procedure,  but  the  appelation  "Newton's  method" 
came  to  be  applied  to  this  whole  type  of  approximations  through 
the  writings  of  Lagrange  and  Fourier. 

3.   raphson's  development  of  newton's  method 

The  first  mathematician  to  enlarge  upon  the  new  method  of 
approximating  roots  by  substitution  in  equations  of  infinite  series 
was  Joseph  Raphson.  We  may  even  say  that  he  was  the  first  to 
expand  it  into  a  system.  It  was  he  who  made  Newton's  method 
a  component  part  of  modern  algebra;  for  with  Newton  it  was  an 
accessory  to  the  calculus,  and  was  never  taken  up  in  his  Universal 
Arithmetic.  In  1690  Raphson  published  a  tract  of  forty-eight  pages 
entitled  Analysis  Aequationum  Universalis  seu  Ad  Aequationes 
Algebraicas  Resolvendas  Methodus  Generalis,  &  Expedita,  Ex  nova 
Infinitarum  Serierum  Methodo.  In  1697  he  published  a  second  edi- 
tion with  an  appendix  discussing  the  recently  published  methods 
of   DeLagny   and    Halley.     Though   quoting   from   Wallis  7   as   to 

6  Letter  to  Mr.  Oldenburg,  Oct.  24,  1676,  given  in  Stewart's  commentary  on  Analy- 
sis by  Equations,  p.  358. 

7  Raphson.  Analysis  Aequationum  Universalis,  London,  1697,  p.  50. 
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Newton's  work  in  infinite  series,  he  does  not  even  refer  to  Newton's 
method  of  approximation. 

Raphson's  method  is  eclectic.  He  sets  up  the  example  like  Vieta, 
employs  Harriot's  notation,  and  divides  the  coefficients  and  resol- 
vends  into  periods,  "modo,  quern  ducere  Vieta,  nostratisque  Har- 
riotus  &  Oughtredus."  Like  Oughtred  he  suggests  the  use  of  loga- 
rithms in  finding  the  first  digits  of  the  root;  he  adds  that  it  may  also 
be  done  by  constructing  parabolas  and  by  examining  the  limits  of 
the  roots.  But  the  main  contributing  elements  in  his  method  are 
Harriot's  use  of  canonical  forms  and  Newton's  method  of  substi- 
tution. Employing  and  enlarging  upon  these  ideas  Raphson  builds 
a  canonical  form  for  the  supplement  of  the  root  in  each  of  sixty- 
four  canonical  equations;  so  that  instead  of  substituting  anew 
in  each  equation,  as  in  Newton's  example,  he  could  substitute  the 
successive  approximations  in  one  canonical  form.  The  "Newtonian 
divisor"  in  these  forms  is  a  misnomer  and  should  be  credited  to 
Raphson.  The  name  given  to  the  collection  of  forms,  "Canones 
Directorii,"  is  the  same  as  that  used  by  Harriot. 

Only  the  two  equations  a3  =  d  and  ba  —  a3  =  c,  representing 
pure  and  affected  powers  respectively,  have  the  canonical  forms 
for  their  supplements  derived  in  full.  Their  derivations  are  as 
follows: 8 

(i)  If  a3  =  d,  and  a  =  g  +  x,  then  we  have  g3  +  $g2x  + 
3gx2  +  x3  =  d.    Rejecting  the  powers  of  x  above  the  first,  as  being 

d  —  g3 

comparatively  small,  we  have  the  supplement  x  — — . 

322 

(2)  If  ba  —  a3  =  c,  and  a  =  g  +  x,  then  "e  Theoremata  Vietae," 
bg  —  g3  _j_  (ft  _  3g2)T  _  3gx2  _  xz  _  Cm    Again,  rejecting  the  powers 

c  _j_  gi  _  0gm 

of  x  above  the  first,  we  have  tor  the  supplement  x  =  — ■ — 

b  -  3g2 

Then  he  proceeds  to  prove  the  convergency  of  the  series. 

We  append   the  following  problem  9  solved   by   this  last   form : 

"Numeris  77284  a  -  aaa  =801   312  8 

8  Raphson,  pp.  5-6,  8. 

9  Raphson,  p.  19. 
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-ru „  _c  +  ggg-bg 


b  -  3gg 

g  =  200 

c  +  ggg  -  bg  = 
+  b  -  2>gg  = 

+  626328 
—  42716)  +  626328(—  14  =  x 

200 

-14 
g   =   186 

£  +  ggg  ~  bg  = 
+  b  -  3,gg      = 

+  143  160 
-  26504)  +  143  i6o(-  5.4  =  x 

186 

-54 

g  =   180.6 

c  +  ggg  -  bg  = 
+  b  -  3gg  = 

+  16152.216 
-  20565.08)  +  i6i52.2i6(-  .7854  =  x 

180.6 

-7845 

g  =  179.8146  vel  179.814  =  g 
c  +  ggg  -  bg  =  +  322.299  405  144 

+  b  -  sgg  =       -  19715.224)   +  322. 299405 1 44 (-  .016348  =  x 

179.814 
—  .016  348 
a  =  179.797652." 

The  tract  contains  the  solutions  of  thirty-four  numerical  problems 
taken  mostly  from  geometry  and  trigonometry.  One  problem,  in- 
deed, deals  with  interest  and  is  undoubtedly  the  first  instance  of 
Newton's  method  being  used  in  the  solution  of  exponential  equa- 

0   —    £365 

tions.10     Using  x  = —  as  the  canon  form  for  the  supplement 

365  g  364 
of  the  equation  a365  =  b,  he  solves  quite  accurately  the  equation 
a365  =  $1.06. 

4.     THE  APPROXIMATIONS  OF  DELAGNY 

Thomas  Fantet  DeLagny  (1 660-1 734)  spent  over  forty  years  in 
work  on  the  approximation  of  roots.    Observing  that  the  expression 

10  Raphson,  p.  36,  problem  29. 
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=  \Za2  +  b  has  a  rational  approximation  z  =  a  +  —  and  a  radi- 

2a 


2—  y  4  — 


cal  approximation  z  =  -  +  4/  —  +  6,  he  sought  to  obtain  similar 


expressions  for  z  =  \/a3  +  6.  These  he  found  to  be  -  +  4/  — h  —  > 

2        r     4       3a 


2  r       A 


.    b  —  I        ,               a6  +  a  a& 

2  >  -  +  4/  -  +  and  a  -\ >  z  >  a  + 


2         y    4  3a  3a2  +  b  +  1  3a2  +  6 

These  results  were  published  in  1691  12  and  in  1697  he  incorporated 
them  in  the  third  part  of  his  Nouveaux  Elemens  d'Arithmetique  et 
d'Algebre,  and  later  in  his  Analyse  Generate  (1733).  He  derived 
similar  expressions  for  higher  powers.  The  process  itself  we  shall 
discuss  in  connection  with  Halley's  general  formula.  But  the 
reader  will  be  interested  to  know  that  in  finding  the  roots  of  large 
numbers,    DeLagny    preferred    to    use    the    formula    \/a3  +  b  = 

a  +  rather  than  the  Hindu  method.   Thus  he  finds  the  cube 

3aa  +  b 


root    of    817,400,375    as    follows:    ^817  =  ^729  +  88  =  9  -^-; 

multiply  this  by  100,  since  the  desired  root  has  three  digits;  the 
result  is  935. 

DeLagny's  method  of  attack  in  solving  both  pure  and  affected 
equations  may  be  seen  from  his  solution  of  the  equation  x3  =  a?x  +  b.u 
Let  x  =  a  +  y-  Then  x3  =  a3  +  3a2y  +  3ay2  +  y3;  also  x3  = 
a3  +  a2y  +  b.    Therefore  b  =  i,ay°  -f-  2a2y  +  y3-    Neglecting  y3  as 

being  very  small,  3ay2  +  2o?y  —  b  =  o,  or  y  = h  4/  -  +  — • 

3         v    9       3a 


Consequently,  x  =  —  +4/  —  +  —  •    Replacing  b  by  b  —  y3  in  the 
3         y     9        3a 

11  The  incorrect  data  given  in  Cantor,  III  (1898),  p.  115,  seem  to  be  taken  from 
Halley's  article  in  the  Philosophical  Transactions  for  1694. 

'2  DeLagny,  Nouveaux  Clemens,  Paris,  1697,  pp.  303-6.  Analyse  Generate,  Paris, 
1733.  PP-  036,  342-46,  362-94.    Journal  des  Scavans,  May  14,  1691,  pp.  297-98. 

13  DeLagny,  Nouveaux  itlemens,  pp.  473-74. 
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last  expression  gives  x  =  —  -f  A/ \-  .    The  formulas  for  x 

3  y     9  3a 

and  y  are  canonical  expressions  for  the  successive  roots  and  supple- 
ments of  the  equation  x3  =  a?x  +  b.    Raphson's  canonical  form  in 

C  ~\~  02  -\~  g3 

the    equivalent    equation    a3  —  ba  =  c    is    a  =  g  + 


3£2  -  b 


different  in  form  and  derived  by  a  different  process.  In  DeLagny's 
illustration,  x3  =  7569*  -f-  240,903;  then  a2  =  7569,  a  =  87,  b  = 
240,903;  this  gives  x  =  58  +  V1764  =  100  for  a  first  approxima- 
tion, x  =  58+4/  1755 — -  for  a  second  approximation,  and  so  on. 
v  261 

DeLagny's  "Method  of  Mediation"  resembled  the  more  effective 
Rule  of  Mean  Numbers  given  by  Chuquet;  except  that  DeLagnv 
used  the  arithmetical  mean.14 

Two  points  about  DeLagny's  work  are  worthy  of  special  men- 
tion. He  stood  alone  among  his  contemporaries  in  his  insistence  on 
a  system  of  gauging  the  error  of  the  approximation.  "Every  method 
is  useless,"  he  says,  "if  it  be  not  accompanied  by  another  method 
which  gives  the  limits  of  error,  be  it  par  exces  or  par  defaut."  15  And 
it  was  his  pride  that  in  the  Triangle  du  Rapport  he  had  an  unfailing 
scheme  for  detecting  the  limits  of  error. 

The  other  point  relates  to  his  ambitious  scheme  of  facilitating 
the  solution  of  higher  degree  equations  by  means  of  tables.  For  this 
purpose  he  worked  out  three  elaborate  sets:  (1)  tables  for  solving 
numerical  equations  directly,  covering  forty-nine  large  pages;  (2)  the 
"Trapese  Logarithmique,"  16  for  facilitating  by  means  of  arithmetical 
progression  the  solution  of  equations;  (3)  the  "Triangle  du  Rapport" 
for  expanding  the  roots  into  series.  These  tables  together  with 
Raphson's  Canones  Directorii  might  have  come  to  fill  a  practical 
need,  as  do  our  integral  and  logarithm  tables  of  to-day,  had  they 
not  been  made  obsolete  by  the  discoveries  of  Halley  and  Taylor. 

DeLagny's  style  was  burdened  with  much  verbiage  and  his 
classifications    and    arrangement    lacked    sharp    lines.     That   may 

14  Nouveaux  Elemens,  pp.  507-11. 

15  Analyse  Generate,  p.  546. 

"  Analyse  Generate,  pp.  320-323;  see  also  Memoires  dr  I' Academic  for  1  705  and  1706. 
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account  somewhat  for  the  comparatively  slight  influence  of  his 
voluminous  work,  some  parts  of  which  are  very  scholarly.  Also,  if 
he  had  lived  in  England  instead  of  in  France,  his  influence  might 
have  been  more  widespread.  At  any  rate,  the  forty-eight  pages  of 
Raphson  left  a  greater  immediate  impress  on  approximation 
methods  than  did  the  twelve  hundred  pages  of  DeLagny. 

5.    halley's  invention  of  general  FORMULAS 

Dr.  Edmund  Halley  (1656-1742),  the  English  astronomer,  had 
learned,  through  a  friend,  of  DeLagny's  results  for  v/a'i  +  ^and 
v  ah  +  b,  and  this  inspired  him  to  investigate  approximation  meth- 
ods for  such  expressions.  It  is  amusing  to  know  that  a  garbled 
version  of  these  results  had  reached  Halley.    For  the  cubic,  for  in- 

a    1      .  /a1       h-i 


stance,    DeLagny   found    that    -  +   A  /  _  _|_ .    <  \/  a3  _)_   0  < 


2 


4  3« 


2         V    4        la  3a3  +  o 


4        sa  3a3  +  b  3a3  +  b  +  1 


a 


/■: 


It  reached  Hallev  in  the  form  a   +    <    Va3  4-  b   <  ~  + 

3a3  +  b  2 

_)__..       As    no   demonstration    accompanied    the    report    of 

DeLagny's  work,  he  derived  these  statements  by  a  method  of  his 
own;  this  was  analogous  to  that  of  the  French  scholar  and  became 
the  basis  of  his  later  investigations.17  He  published  his  results  in  the 
Philosophical  Transactions  for  1694  in  an  article  entitled  "A  ne>v 
exact,  and  easy  method  of  finding  the  roots  of  any  equations  gener- 
ally, and  that  without  any  previous  reduction."  In  this  article  he 
differs  from  Raphson  and  Halley  in  that  he  gives  frank  and  generous 
recognition  to  the  work  of  his  contemporaries. 

His  approximation  for  V a3  +  b  is  as  follows: 

(1)  Suppose  a3  +  b  to  be  a  non-cubic  number;  also  suppose  that 
a3  -f  b  =  (a  +  e)8,  approximately,  so  that  a*  <  (a  +  e)3  <  (a  +  i)3. 

Then  (a  -f-  e)3  =  a3  +  y&e  +  2>ae-  +  e3  =  a3  +  b  and  (a  +  i)3  = 
a3  -+-  3a-  +  yi  +  1  ;  hence  c:i<  1. 

Philosophical  Transaction  of  the  Royal  Socirty  <>./  London.  \<>.  210  (1694),  p.  136; 
abridged  English  translation  by  Mutton.  Shaw,  and  Pearson,  III  (1683-1694),  pp.  640- 
Y).    Hereafter  referred  to  as  Phil.  Trans.  (Abr.). 
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Also   b  =  3a2e  +  3ae2  +  es;   or   b  =  i,a2e  +  3ae2,   approximately. 
Since  a?e  is  so  much  greater  than  ae2,  he  thinks  of  e  as  equal  to 
b_ 

3^' 

h  h 

From  b  =  3a2 ]e  +  3a?2,  he  gets  e 


3a3  +  3ae 
a& 


3aS  +  3a  *>3a 


(-) 

Ma2/ 


3a3 +  6 


a& 


Therefore  Va3  +b  =  a  +  e  =  a  + —  •    This  is  the  rational 

3a3  +  b 

form,  which  by  his  handling  of  e  is  slightly  less  than  a  +  e. 

(2)     Since     b  =  2,a2e  +  2>ae3,     it     follows     that     e2  -\-  ae  =  — ; 

completing   the  square  and   finding  roots  gives   the   radical   form 


a  +  e  =   -  -\-  A/  —  -\-  — ,  which  is  a  little  too  large.  Similar  proc- 
2         f    4         3a 

ab  ^       / ^   #    ,         /a2         h 

esses  give  a <   Va3  -  6  <  -  +  \     -  -  — .   The  other 

results  reported  by  DeLagny  are  then  verified. 

So  far  Halley  has  only  re-discovered  DeLagny's  results.  But  now 
he  takes  the  important  step  which  constitutes  his  own  real  contribu- 
tion to  facilitating  the  operation  of  Newton's  method.  First  he 
derives  a  general  form  for  all  pure  powders;  thereupon  he  extends  the 
method  of  computation  as  well  as  its  generalization  to  affected  equa- 
tions. It  was  this  important  feature  that  led  to  Taylor's  important 
discovery.  Halley  says:  "While  I  was  considering  this  matter,  I  fell 
upon  a  general  method  of  forms  for  any  power  whatever,  which 
is  elegant  enough,  and  which  I  cannot  prevail  on  myself  to  conceal, 
since  in  the  higher  powers  they  at  least  triple  the  known  figures  of 
the  root.  And  these  forms,  both  rational  and  irrational,  proceed  as 
follows : 
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\/a3  +  b 


0  ,        /i               b  ,  ab 
-a  +  A/  -aa  -\ 5,  or  a  + 

1  K     1  ia 

=  -  a  +  A/  -aa  H ,ora  + 

2  T    4  3« 


v/a2  +  6  =  -c+/l//-aaH 5,  or  a  + 

1  V     l  la  2a2  +  -b 

ab 


3<*3  +  f* 


2  .  /i  6  ,         ab 


1/ 


vV  +  b  =  -a+A/-aa  +  — ,,  or  a  + 


3  ^    9  6a2  4a4  +  3_/; 


!»+i/^ 


b  ab 


</ab  +  b  =  -a+A/—  aa-\ -,  or  a  + 


5«5  +  |i 


°VI 


^a8  +  6  =  -  a  +  4/—  aa  H ,,  or  a  + 


5  V    25     '        '5«<  6flt  +  56 

2 

And  so  on  for  the  higher  powers." 

For  the  purpose  of  extending  this  method  to  affected  equations  of 
any  degree  he  builds  up  what  he  calls  a  "general  analytical  spe- 
culum."   This  is  a  table  of  powers  to  aid  in  the  transforming  of  the 

equation     cz  +  dz2  +  /z3  + =0    by     letting    z  =  a  -f  e. 

The  resulting  coefficients  of  e,  e1,  e3,  etc.,  are  represented  by  s,  t,  u, 
etc.    Following  is  his 

Table  of  Powers 
s  t  u  w  x  y 

cz    =  ca    -f  ce 
dz2  m  da2  +  2dae    +    de2 

fz*  =  /a3  +  ?>fa2e  +  ifae2  +  //e3 

gz4  =  gaA  -f  4ga3e  -f  6gaV  -f  4gae3      +  ge4 

/zz5  =  ha>  +  5^a4e  -f  ioha3e2  +  10/zaV  -f  5/me4      +  heb 

kz*  =  fca6  +  6ka*e  +  isifcaV  +  20ka3e3  +  1  skate*  +  5fcae6  +  ke6 

The  modern  reader  recognizes  the  sums  of  the  coefficients  in  the 
different  columns  to  be  the  same  as  in  the  expansion /(z)  =  f(a  +  e) 


f(a)   {-f'(a)e  + 


f^{a)  g2  +  f'"(a)  ea  +  ....  But  Halley  denotes  the 


2!  3! 
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result  in  the  form  analogous  to  that  which  he  uses  for  pure  powers, 

viz.,  b  ~  se  +  te2  +  ue3  +  wtA  -f  xeh  -f  ye*  +  ,  where  6  =  /(s)  —/(a). 

Operating  with  this  equation  he  finds  that  the  approximate  roots 

of  s  are:  s  =  a  ±  — — — ,    in    the    rational    form,    and   z    =    a    ^ 
s-  ±  tb 


\/\  <'*l» 


I 

s± 

in  the  radical  form. 


Among  his  numerical  illustrations  is  the  equation  z4  —  80s3  + 
1998s2  —  14937s  +  5000  =  o,18  taken  from  Wallis.  Raphson  uses 
the  same  illustration.  Even  a  cursory  examination  shows  the  merits 
of  Halley's  General  Analytical  Speculum.  The  arrangement  is  still 
used  in  some  modern  texts  on  the  theory  of  equations. 


6.   taylor's  application  of  the  calculus  to 
halley's  coefficients 

The  next  improvement  came  from  the  field  of  the  differential 
calculus.  Brooks  Taylor  (1685— 1731)  had  shown  in  his  Methodus 
Incrementorum  that 

f(a  +  e)  =  f(a)  +  f  (a)  e  +  f— ^  e?  +^— ^  e'+ (Taylor's 

2!  3! 

Series). 

In  an  article  in  the  Philosophical  Transactions  for  171719  entitled 
"An  attempt  towards  the  improvement  of  the  method  of  approxi- 
mating, in  the  extraction  of  the  roots  of  equations  in  numbers,"  he 
shows  the  identity  of  the  coefficients  in  his  series  and  the  coefficients 
in  Halley's  General  Analytical  Speculum.  Since  his  series  is  general, 
applicable  to  polynomials  as  well  as  binomials,  to  fractional  as  well 
as  integral  exponents,  to  transcendental  as  well  as  algebraic  expres- 
sions, he  observes  that  Halley's  work  can  easily  be  extended  into 
these  fields  by  aid  of  the  said  series.  He  illustrates  this  by  solving 
the  equations  ( v'2  +  0'1  +  y  —  16  =  O  and  log  y  =  0.29. 

18  Phil.  Trans.  (Abr.)i  P-  648;  Wallis,  Algebra,  pp.  233-35;  Raphson,  pp.  25-26. 

19  Phil.  Trans.  (Abr.),  VI  (1713-1723)  pp.  299-304;  Phil.  Trans.,  Vol.  30.  pp.  610- 
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7.     FURTHER  SIMPLIFICATIONS  OF  NEWTON'S  METHOD 

In  modern  text-books  we  find  the  Newtonian  method  of  approxi- 
mation stated  thus:  if/(.r)  =  o,  and  a  is  a  first  approximation,  then 

x  =  a  —  — — -  is  a  second  approximation.    Newton  and   Raphson 
f'(a) 

had  used  these  values,  but  obtained  them  by  transforming  the  equa- 
tion. Halley  used  the  quadratic  part  of  the  transformed  equation; 
Newton  did  the  same  at  times.  Brooks  Taylor  also  used  the  quad- 
ratic,  but   obtained   his  coefficient    by   differentiation.     The    first 

writer  to  use  the  linear  differential  formula  x  =  a  —  was  Thorn- 

f'(a) 

as  Simpson  in  his  essay,  "A  new  Method  for  the  Solution  of  Equa- 
tions in  Numbers"  (1740).20  He  states  the  formula  rhetorically  and 
illustrates  it  by  this  example:  f(x)  =  30OX  —  x3  —  1000  =  o;  x  = 
a  +  A;  by  trial  a  is  found  to  lie  between  3  and  4;  try  a  =  3,5  ]f(x) 
=  300  -  3.T2;  /(a)  ^/'(a)  =  A;  f^.5)  =  7.125;  f'(2>-5)  =  263.25; 
A  =  0.027;  3-5  -  0.027  =  3473- 

Proceeding  as  before,  letting  a  =  3.473,  we  have  x  =  3.47296351. 

Simpson  uses  Newton's  method  in  solving  simultaneous  equations. 
His  fifth  example  is  a  set  of  exponential  equations:  x*  -f-  y1  —  1000 
=  o  and  xy  -f-  yx  —  100  =  o. 

De  Courtivron  (1744),  Euler  (1755),  and  Lambert  (r.  1770)  also 
applied  the  calculus  to  the  solution  of  higher  equations.21  In  his 
algebra  w  (1770),  Euler  discards  the  higher  powers  in  both  dividend 
and  divisor  during  the  substitution  instead  of  carrying  them  along 
until  the  division  takes  place.  Thus  in  solving  the  equation  .r:i  + 
ax2  +  bx  +  c  =  o,  he  sets  x  —  n  —  p,  x2  =  ri-  —  inp,  x3  =  n*  — 

•>,.    r         ,,         c  .        ,     c    ,     ,        »3  +  an2  +  bn  -f  c 

,VJ-p;  from  these  lour  expressions  he  finds  p  = • 

3«2  +  2an  +  b 

Edward   Waring    (1770)    extended    the    Newtonian   approximation 

process  into  the  realm  of  complex  roots.    The  functional  expression 

x  =  a  —  — —  is  due  to  Lagrange23  (179H). 
/  (a) 

-''  Simpson.  Essays  on  Mathematics,  London,  1740,  p.  81. 

a  Euler.  Instilutiones  calculi  differenlialis,  sections  224.  234,  235;  Histoire  de 
I' Academic.  1744.  p.  405—14. 

-  Eukr,  Algebra,  I  (1770),  Chap.  16. 
•  Lagrange,  Traite  de  la  Resolution  ties  Equations  Nutneriques  |  179S).  Nop-  X  I. 
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i.  The  method  of  solving  equations  of  higher  degree  by  suc- 
cessive substitutions  in  dcived  equations  was  originated  by  Newton 
(1669)  incidental  to  his  getting  integral  expressions  for  his  work  in 
areas.  He  performed  a  transformation  for  every  new  supplement  to 
the  root.  His  method  was  first  made  public  in  the  algebra  of  Wallis 
(1685). 

2.  Raphson  (1690)  worked  Newton's  method  into  a  system. 
Following  the  example  of  Harriot  he  derived  canonical  forms  to 
simplify  calculation.  In  his  plan  only  one  transformation  and  one 
divisor  are  needed  for  each  equation. 

3.  DeLagny  (1691)  derived  notable  approximations  in  the  solu- 
tion of  pure  equations  by  a  special  adaptation  of  Newton's  method, 
and  worked  out  an  elaborate  set  of  tables  for  the  solution  of  numer- 
ical equations. 

4.  Halley  (1694)  generalized  DeLagny's  method  for  pure  equa- 
tions and  extended  it  to  affected  equations;  in  his  table  of  powers  he 
gave  a  general  formula  for  the  coefficients  of  transformed  equations. 
This  made  obsolete  Raphson's  canonical  forms  and  DeLagny's 
tables. 

5.  Taylor  (1717)  discovered  that  Halley's  coefficients  were  the 
same  as  the  coefficients  in  Taylor's  series.  This  made  possible  the 
direct  use  of  the  calculus  in  finding  roots. 

6.  Our  present  standard  linear  approximation  formula  x  =  a  — 

was  first  used  by  Simpson  (1740);  it  received  its  conventional 

form  through  the  writings  of  Lagrange  (1798). 


VIII 

CERTAIN  EFFECTIVE  BUT  NON- 
PRACTICABLE  METHODS 

I.     THE  ASCENDANCY  OF  NEWTON'S  METHOD 

Newton's  method  of  approximating  roots  by  substituting  values 
in  equations  of  an  infinite  number  of  terms  did  not  immediately 
upon  its  publication  in  1685  supersede  Vieta's  method  of  finding 
roots  by  evolving  the  successive  digits.  The  latter  was  still  used  by 
Dechales  in  1690  and  appears  as  late  as  1702  in  a  new  edition  of 
Oughtred's  Clavis.  It  seems  as  if  the  transition  from  the  one  to  the 
other  must  have  been  facilitated  by  the  eclectic  character  of  Raph- 
son's  arrangement.  Halley's  general  formula  increased  the  popular- 
ity of  the  new  method.  Its  ascendancy  is  clearly  shown  in  an  article 
written  by  Wallis  in  1694.  After  Taylor's  discovery  (1717)  it  had  no 
serious  rival  old  or  new,  until  the  publication  of  Horner's  method  in 
1819.  But  new  methods  sprang  up  to  challenge  it  in  this  period,  and 
we  shall  describe  the  three  that  gave  the  most  promise. 

2.     ROLLE'S  METHOD  OF  CASCADES 

A  considerable  contribution  to  the  location  of  roots  was  made  by 
Michel  Rolle  in  his  two  works  Traite  d'Algebre  (1690)  and  Demon s- 
tratio  d'une  Methode  pour  Resoudre  les  Egalitez  de  toutes  Degrez 
(1691).  He  demonstrated  that  there  cannot  be  more  than  one  real 
root  of  /(.%•)  =  o  between  two  successive  real  roots  of  f{x)  —  o 
(Rolle's  Theorem).  He  invented  a  method  of  finding  roots,  called 
the  "Method  of  Cascades";  it  consisted  in  locating  the  roots  by  first 
finding  the  limits  ("hypotheses")  of  the  roots  of  the  different  cas- 
cades. Between  the  roots  of  each  cascade  lie  the  roots  of  the  next 
higher  cascade,  and  so  on.1  It  was  an  anticipation  of  the  method  of 
using  successive  derivatives  (Rolle's  cascades)  now  used  in  calculus 
for  locating  roots. 

With  Rolle  an  equation  is  "piepared"  when  it  has  the  form 
x"  —  ai.r""1  +  a2xn~2  —  azxn~A  + =0  where  the  a's  are  posi- 

1  Traite  d'Algebre.  p.  124.  p.  103;  p.  1 18;  Demonstratio,  article  VI  (New  York  Public 
Library);  Cajori  in  Bibl.  Math.  XI  (1910-11).  p.  300  313. 
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live  integers.    Then  the  upper  limit  of  the  root  is  found  by  dividing 

the  largest  negative  coefficient  by  the  highest  exponent,  and  adding 

unity  and  as  much  more  as  necessary  to  make  the  result  an  integer. 

His  Method   of  Cascades  he   illustrates   thus:2 

Let  x3  — 57-X"2  -f-  936*  —  3780  =  o  (third  cascade). 

Multiply  the  terms  of  the  cascade  by  parts  of  the  piogression  3,  2, 

1,  o,  respectively;  divide  by  x.  Then  3X2  —  114.%  -f-  936  =  o  (second 

cascade). 

Multiply   by  2,   1,  o;   divide  by   2x.   Then    3.r    —   57    =    o    (first 

cascade) . 

The  "hypotheses"    for    the    second    cascade  are  o,    19  and   39   (or 

H-4  +  1). 
3 

The  "hypotheses"  for  the  first  cascade  are  o  and   iq. 

Between  o  and  19,  the  second  cascade  is  found  to  have  12  as  a  root; 

and  another  root,  21,  between   19  and  39. 

Finally  for  the  third  cascade  the  extreme  "hypotheses"  are  o  and  3781 


or 


3780 


+  1  j  ;  o,    12,   26,   and   3781    are   its  four   "hypotheses. 


The  roots  of  the  third  cascade  (that  is,  the  given  equation)  can 
then  be  found  to  be:  6  between  o  and  12;  21  between  12  and  26;  30 
between  26  and  3781. 

3.    THE  METHOD  OF  RECURRING  SERIES 

De  Moivre  (1667-1754)  first  conceived  of  series  in  which  each  co- 
efficient bears  a  given  relation  to  certain  of  the  preceding  coeffi- 
cients; such  series  he  called  lecurring  series  (1720).3  Daniel  Ber- 
noulli (1700-1782)  first  used  them  for  approximating  roots  (1728). 
Following  is  his  process:4 

Take  such  an  equation  as  1  =  ax  +  bx2  +  cxs  -f-  dx4.  Select 
arbitrarily  the  four  terms  A,  B,  C,  D.  Then  E,  F,  G,  H,  etc.,  con- 
stitute a  recurring  series  if 

2  Traite  d'Algebre,  pp.  127-28. 

3  Phil.  Trans.,  Vol.  32  (1722),  pp.  162-78;  DeMoivre,  Miscellanea  Analytica, 
London,  1730,  Bk.  I,  II. 

4  Cantor.  Ill  (1898),  pp.  621-22. 
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E  =  aD  +  bC  +  cB  +  <IA, 

F  =  aE  +  bD  +  cC  +  dB, 

G  =  aF  +  bE  +  cD  +  d( ', 

H  =  a-G  +  bF  +  cE  +  aD; 

E  F 
and   the  approximate  values  of  x  are-,-,  etc. 

F  G 

In  solving  the  equation  1  =  —  2.Y  +  5.Y1  —  4.x3  +  .\\  Bernoulli  sets 
.4  =  B  =  C  =  D  =  1,  and  obtains  the  series 
0,2,  -7,25,  -93.341.  -1254. 

This   gives   -,  —  =,  —    — , - —    as    the    increasingly   close 

2         7         25  1254 

approximations  of  the  root.    The  last  value,  when  substituted    in 
the  given  equation  gives  1  =  0.999487,  a  very  small  error. 

Bernoulli  gave  no  analytical  proof  for  the  validity  of  his  process. 
That  was  supplied  by  Euler  (1748)  and  later  by  Lagrange  (1798).* 

A  considerable  improvement  was  made  by  Euler  in  his  algebra 
(1770).  He  shows  the  general  relation  that  must  exist  for  the  dif- 
ferent powers,  and  thus  immediately  gets  a  recursion  formula  for 
any  equation  under  consideration.6  The  recursion  formula  performs 
a  function  similar  to  that  of  Raphson's  canonical  form  or  of  Halley's 
general  formula. 

n        y 

Suppose  our  series  p,  q,  r,  S,  t,  etc.,  is  such  that  -,  -,  etc.,  express 

P    (J 

approximate  values  of  x.   Then:      =  x\-  =  x\  multiplving,  -  =  .v'J; 

p  q  p 

s  s  .     .  t 

also,     =  .v;  hence-  =  .v3;  similarly  -  =  .r1,  etc.     Not  only  do  we 
r  p  p 

have  a  substitution  table  for  powers,  but  we  have  one  with  the  first 
term  as  the  common  denominator. 

As  an  illustration,  let  us  substitute  these  values  in  the  equation 

x3  =  Xs  4-  2.v  4-  1.    This  gives  s  =  r  -\-  2q  4-  p-    Using  this  as  a  re- 

6  Euler's   InlroJuctio   in    Analysin    1  nfuiilorum,    (hap.    17;   Lagrange,    Traiti   Jf  la 
Resolution  des  Equations  Numeriques,  Note  VI. 
*  Euler's  Algebra,  Vol.  I.  Chap.  16. 
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cursion  formula  for  5  in  terms  of  r,  q,  and  p,  we  can  arbitrarily  select 
these  to  be  0,0,  1.   We  then  get  as  our  series  0,0,  1,  1,  3,  6,  13,  28 

60,  129,  etc.    Testing  for  x  =  —  in  the  original  equation  we  get  an 
28 

j  -j  i  I2Q 

error  of  — .    A  closer  result  would  be  given  by  x  =  — -. 
343  60 

Not  all  equations  admit  of  the  application  of  the  method  of 
recurring  series.  If  no  common  limit  results,  the  equation  must  first 
be  transformed.  This  happens  particularly  if  the  second  term  is 
lacking.  Thus  foi  the  equation  x2  =  2  we  get  the  recurring  series 
1,  1,  2,  2,  4,  4,  8,  8,  etc.  A  similar  situation  arises  for  x3  =  2.  To 
obviate  this,  we  let  x  =  y  —  1.  In  the  resulting  equation  y3  = 
2>J2  ~  3y  +  l  we  get  the  recursion  formula  5  =  2>r  ~  33  +  P>  and 
this  gives  us  the  series  0,0,  1,  3,  6,  12,  27,  63,  144,  324.   This  gives  us 

y  =  —  and  x  =  -.    Checking  will  reveal  an  error  of  —  . 
144  4  64 

4.     THE  METHOD  OF  CONTINUED  FRACTIONS 

Newton's  method  had  one  serious  defect:  it  was  inoperative  if 
two  roots  were  close  together,  for  then  the  series  would  not  converge, 
and  might  after  a  while  actually  diverge.  Besides  this,  Lagrange 
pointed  out  two  other  defects:  the  uncertainty  as  to  the  exactness 
of  each  fresh  correction,  and  the  failure  to  give  a  commensurable 
root  in  finite  terms.  Lagrange  announced  (1767)  a  new  method 
which  was  free  from  these  defects.7  It  consisted  of  three  parts:  (1) 
a  scientific  method  of  finding  the  integral  part  of  the  root  instead  of 
the  empirical  methods  hitherto  employed;  (2)  a  rule  for  separating 
the  roots;  and  (3)  a  technique  of  approximation  by  the  use  of  con- 
tinued fractions.    This  last  may  be  explained  as  follows: 

Suppose  that  in  the  equation /(x)  =  o  we  have  found  that  there  is 
only  one  root  between  the  integers  p  and  p  -f-  1 ;  for  convenience  let 

it  be  real  and  positive.   Then  x  =  p  +  -;   and   y   is   greater   than 

y 

unity.   Substituting  the  value  of  x  \nf{x)  gives  an  equation  F{y)  =  o. 

'  Lagrange.  Traile,  Chap.  I,  III,  IV  and  Notes  IV,  V. 
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Since  y    >    1,  find   its  integral  value  by   the  prescribed  method. 

Suppose  it  lies  between  the  integers  q  and  q  +  1.    Then  y  =  q  -\-  -  ; 

z 

and  so  on.    Then  x  =  p  +  - 


2  +     ] 


r  +  .... 

Lagrange  illustrates  his  method  by  numerical  examples.  One  of 
these  is  Newton's  equation  x3  —  2*  —  5  =  o.8  By  his  methods  of 
limits  and  separation  he  finds  one  root  between  2  and  3.    Let  x  = 

2  -f-  -;   this  gives   the  equation  yz  —  ioy2  —  6y  —  1  =  o.      Again 

y 

applying  the  method  of  limits  and  separation  he  finds  that  y  lies 
between  10  and  II.  Similarly  for  y  =  10  -+-  -  he  derives  the  equa- 
tion 6123  —  94s2  —  20z  —1=0;  here  I  <  s  <  2. 

Letting  2=1+-     gives     5411s  -f  25W2  —  89W  —  61  =0:    here 
u 

1   <  u  <  2. 

A  continuation  of  this  process  gives  the  series  2,  10,  1,  1,  2,  I,  3, 
I,  1,  12. 

Hence  x  =  2  +  I 


10  +  1 


1  +  1 


1  +J 

2  + 


Evaluating  this  expression  by  known  methods  gives  the  approxi- 

2    21    23    44    in     115    576    731    1307    16415 

mating  series:  -,—,—,—,  ,  — -,  ^—,  1^~,  -^-L,  — 3-£,  etc. 

1    10    11    21     53      75     275    349     624      7837 

»  Traitl  (1798),  Chap.  IV,  p.  29. 
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That  is,  x   =     -— LJ?  =  2.09455149;    and    the    error    is    less    than 
7837 
1 

0.0000000163. 


(7837) : 


Though  none  of  these  methods  gained  popular  usage,  they  gave 
rise  to  by-products  that  were  incorporated  by  algebraists  into 
Newton's  method.  Among  these  were  Rolle's  theorem  for  locating 
roots  and  the  method  of  Lagrange  for  finding  the  integral  parts  of 
the  root,  as  well  as  his  rule  for  separating  roots.  Waring  had 
already  (1762)  formulated  such  a  rule  but  it  was  very  little  known 
until  Lagrange  published  his  method  of  squared  differences. 

5.     SUMMARY 

1.  Between  1685  and  1819,  three  potential  rivals  of  Newton's 
method  sprang  into  being:  Rolle's  Method  of  Cascades  (1690), 
Bernoulli's  (1728)  and  Euler's  (1748,  1770)  method  of  recurring 
series,  and  Lagrange's  method  of  continued  fractions  (1767). 

2.  None  of  these  methods,  though  scientific  and  effective,  became 
generally  adopted,  and  Newton's  method  was  supreme  from  about 
1694  to  1819. 


IX 

HORNER'S  METHOD;  SIMILAR  METHODS  BY  RUFFINI 

AND  BY  THE  CHINESE  OF  THE 

THIRTEENTH  CENTURY 

I.  A  GENERAL  VIEW 

No  method,  not  even  that  of  Vieta,  had  been  so  universally 
adopted  as  that  of  Newton.  Of  its  two  latest  rivals  the  method  of 
recurring  series,  though  as  presented  in  Euler's  Algebra  not  par- 
ticularly laborious,  never  had  its  possibilities  exhausted;  and  La- 
grange's method  of  continued  fractions,  though  effective  and  perfect 
in  theory,  was  too  laborious  for  popular  adoption.  At  the  beginning 
of  the  19th  century,  therefore,  Newton's  method  held  undisputed 
sway.  Its  one  drawback  had  been  its  insecurity.  But  the  success  of 
Waring  and  Lagrange  in  effecting  methods  for  the  separation  of  the 
roots  had  to  a  large  degree  eliminated  this  defect.  The  discoveries 
of  Fourier  and  Budan,  early  in  the  nineteenth  century,  whereby  the 
number  and  nature  of  the  roots  could  be  detected  for  any  interval, 
gave  Newton's  method  a  much  more  scientific  basis.1  The  mechan- 
ism of  its  operation  had  been  so  greatly  improved  by  Halley,  Taylor, 
and  Simpson,  that  in  this  respect  it  has  to  this  day  remained  without 
a  peer. 

It  has  had  only  one  serious  rival.  It  is  that  form  of  approxima- 
tion by  the  evolution  of  digits  made  known  to  the  Western  world  by 
William  Horner  in  1819.  Like  the  methods  of  Vieta  and  Newton,  it  is 
also  a  process  of  exhaustion.  Horner's  process  resembles  that  of 
ordinary  division  or  of  current  methods  of  finding  square  and  cube 
roots.  That  element  was  not  new,  having  been  the  chief  principle  in 
Vieta's  method.  But  Horner's  way  of  doing  this  by  the  gradual 
diminution  of  the  roots  in  a  succession  of  transformed  equations, 
was  new  to  Europe.  To  the  world  it  was  old,  however;  it  had  been 
known  and  used  in  at  least  three  mathematical  centers  before  [819, 

1  As  early  as  1768  J.  Raym.  Mouraille  had  successfully  formulated  a  plan  for  insur- 
ing the  security  of  Newton's  method;  but  his  work  was  unknown  at  the  time,  and  has 
only  recently  been  noticed.  See  article  by  F.  Cajori  in  Bibl.  Math.  (3)  XI  (191  •>•  PP- 
132-37- 
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— China,  Japan,  and  Italy.    In  this  chapter  we  shall  give  a  short 
survey  of  the  different  forms  this  method  has  assumed. 


2.     THE  METHOD  USED  IN  CHINA  AND  JAPAN 

As  early  as  the  first  half  of  the  7th  century  Wang  Hs'iao  worked 
with  numerical  cubic  equations  in  his  Arithmetic  in  Nine  Sections.2 
He  found  the  roots  by  using  a  method  similar  to  that  of  finding  the 
cube  root  of  numbers.  He  attempted  no  equations  of  degree  higher 
than  the  third. 

Not  till  six  centuries  later  do  we  find  any  work  with  equations 
higher  than  cubics.3  Then  we  find  three  men,  living  almost  at  the 
same  time,  systematizing  a  method  for  the  orderly  evolving  of  the 
digits  of  the  root.  The  first,  Ch'in  Chiu-shao,  in  his  Nine  Sections 
of  Mathematics  (1247),  solved  equations  of  the  6th,  7th,  8th,  and 
even  higher  degrees.  He  worked  on  tables,  with  wooden  rods  to 
represent  his  coefficients,  and  arranged  his  work  in  columns  as  did 
the  Hindus  for  roots  of  numbers.  In  writing  his  equations  with  the 
right  hand  member  equal  to  zero  he  anticipated  Harriot  by  over 
three  hundred  years.  Li  Yeh,  the  second  pioneer  in  this  work,  in 
The  Mirror  of  the  Mensuration  of  Circles  (1259),  explained  the 
process  more  fully;  Ch'in  had  chiefly  given  rules.  The  third  member 
of  the  trio  was  Chu-chieh.  His  form  of  the  process,  as  presented  in 
his  Introduction  to  Mathematical  Studies  (1299),  found  its  way  into 
Japan  in  the  seventeenth  century.  The  names  given  to  this  method 
are  quite  suggestive.  Ch'in  Chiu-shao  (1247)  called  it  "Harmoni- 
ously alternating  evolution";  Yang  Hwuy,  another  writer  of  the 
same  period,  called  it  "Accumulating  involution." 

To  show  the  steps  in  this  process,  invented  in  China  in  the 
thirteenth  century  and  introduced  into  Japan  in  the  seventeenth 
century,  we  shall  give  the  solution  of  a  problem  in  which  the  root  is 
12,  taken  from  the  eighth  book  of  Sato  Moshun's  Tengen  Shinan 
(1698).   The  arrangement 4  is  Occidental.   The  solution  is  as  follows: 

'  Yoshio  Mikami,  The  Development  of  Mathematics  in  China  and  Japan,  New  York 
1914;  pp.  53-56. 

a  D.  E.  Smith  and  Yoshio  Mikami,  .4  History  of  Japanese  Mathematics,  Chicago, 
1914;  pp.  49-56;  Yoshio  Mikami,  pp.  63-89. 

4  Given  by  D.  E.  Smith  and  Y.  Mikami.  The  Chinese  arrangement  is  given  in  full  in 
Smith  and  Mikami's  Hist,  of  Jap.  Math.,  pp.  53-55. 
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Given  xi  +    4*3  -  236X2  -    432*  +     11520  =  o 
1+4      -  236      -    432    +     1 1520 
10  140      —    960    —     13920 


1 

14 

10 

-  96 
240 

-  1392 
1440 

—  2400 

I 

24 

10 

144 
340 

48 

—  2400 
2400 

I 

34 
10 

484 

48 
1 152 

0 

I 

44 

2 

484 
92 

1200 

I      46       576 

Among  Japanese  scholars  to  write  on  the  Chinese  method  6 
of  approximation  were  Sato  Seiko  (1666),  Seki  Kowa  (1 642-1 708), 
Sato  Moshun  (1698),  Sakabe  (1803),  and  Kawai  (1803).6 

3.     THE  METHOD  INVENTED  BY  RUFFINI 

In  1804  a  solution  very  similar  to  the  approximation  process  of 
the  Chinese  was  effected  by  Paolo  Ruffini.7  By  means  of  the  differ- 
ential calculus  he  worked  out  a  theory  for  transforming  one  equation 
into  another  with  roots  diminished  by  a  certain  constant.  In  his 
revisions  of  1807  and  1813  he  used  algebraic  processes  for  the  same 
purpose.  His  arrangement  differs  from  that  of  Horner  in  that  the 
coefficients  of  the  transformed  equation  appear  in  the  extreme  right 
hand  column,  whereas  in  Horner's  arrangement  they  appear  in  a 
diagonal  line.  A  single  solution  will  tend  to  show  the  similarity  of 
Ruffini's  work  to  Horner's  process  and  to  that  of  the  Chinese: 

6  Whether  the  process  used  by  the  Chinese  was  actually  the  same  as  that  used  by 
Horner  is  questioned  by  Professor  Gino  Loria,  of  the  University  of  Genoa,  Italy,  in  an 
article  published  in  the  Bollelino  della  Malhesis  for  April,  1920,  and  abridged  by  Pro- 
fessor R.  B.  McClenon,  of  Grinnell  College,  Grinnell,  Iowa,  in  The  Scientific  Monthly 
for  June,  1921.  In  discussing  Ch'in  Chin-Shao's  solution  of  the  equation  -x*  +  763, 
200X5— 40,  642,  560,  000  =  o  he  says:  "We  still  have  too  meager  information  as  to 
the  details  of  the  work  for  us  to  be  able  to  affirm  confidently  that  Horner's  method  was 
known  to  the  Chinese  in  the  13th  century;  we  can  only  say  that  this  method,  or  one 
practically  identical  with  it,  was  known  at  the  time,  and  we  must  await  further  evidence 
before  affirming  or  denying  the  priority  of  the  Chinese  in  its  discovery." 

'  Smith  and  Mikami,  pp.  53,  86,  116,  212-15. 

7  F.  Cajori,  "Horner's  Method  of  Approximation  Anticipated  by  Ruffini,"  in  Bull. 
American  Math.  Soc,  XVII  (1911),  pp.  409-414. 
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"Given  the  equation  4.x5  —  6xA  -f  3#3  —  5*2  —  4*  +  8  =  o;  re- 
quired, to  transform  it  into  another  in  which  y  =  x  —  6.  Operate 
in  the  manner  here  set  forth : 

4. 


6,  3. 

-  5. 

-4.    8, 

4,  18, 

in, 

661,  3962,  23780 

4. 

42. 

363,  2839,  20996 

4. 

66,  759-  7393 

4.   90,  1299 

4.   114 

4 

and  we  obtain  the  transformation 

4.V5  +  1 14/  +  12993'3  +  73933'2  +  20996V  +  23780  =  o." 

For  Newton's  classic  example  x3  —  2x  —  5  =  o  he  finds  the  root 
to  8  decimal  places,  by  this  method. 

Ruffini's  method  remained  unnoticed  by  his  own  people.  It  was 
the  story  of  Vieta  repeated.  But  there  was  no  Harriot  or  Oughtred 
to  bring  the  message  to  England  this  time.  For  Britain  had  her  own 
discoverer  in  this  field.  It  is  a  coincidence  that  every  successful 
approximation  method  for  numerical  equations  has  been  either  in- 
vented or  sponsored  in  Great  Britain. 

4.   horner's  method 

The  theory  and  mechanism  of  Horner's  method  are  so  well  known 
that  we  need  not  discuss  them  here.8  We  shall  review  only  the  out- 
standing points  in  its  history. 

It  is  safe  to  say  that  neither  Horner  nor  Ruffini  knew  of  the  other's 
work,  nor  was  either  acquainted  with  the  methods  of  the  Orient. 
Horner's  paper  was  read  before  the  Royal  Society  of  London  on  July 
1,  1819,  and  was  published  in  the  Philosophical  Transactions  for  that 
year.9 

In  developing  his  theory,  he  used  Taylor's  Theorem  and  Arbo- 
gast's   derivatives    (Rufnni    had    used    ordinary    derivatives),    and 

8  A  thorough  discussion  of  its  analytic  and  operational  phases  is  to  be  Found  in  J.  R. 
Young's  Theory  ami  Solution  of  Algebraical  Equations  of  Higher  Order,  London.  1843, 
Chapter  XII.  A  vigorous  challenge  to  rival  English  claimants  for  priority  over  Hor- 
ner is  to  be  found  in  De  Morgan*s  article  already  referred  to. 

9  W.  H.  Horner,  "A  new  method  of  solving  numerical  equations  of  all  orders  by 
continuous  approximation,"  in  Phil.  Trans.,  Vol.  109  (1819).  pp.  3o8^335- 
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arranged  his  results  in  a  "General  Synopsis"  in  much  the  same  way 
as  Halley  built  up  his  General  Analytical  Speculum.  He  obtained 
his  transformations  by  using  the  calculus,  but  says  that  algebraic 
processes  can  be  used  for  this  purpose.  Like  Rufnni  he  used  ordi- 
nary algebra  in  his  revisions  of  1830  and  1845. 

His  first  illustration  is  the  classic  example  of  Newton.  The  work 
is  arranged  differently  from  that  which  we  see  in  text-books  on  al- 
gebra. His  schema  resembles  Wallis's  arrangement  of  Vieta's  process 
and  our  present  arrangement  in  finding  the  cube  root;  he  even 
divides  the  resolvend  into  periods. 

Horner's  method  met  with  quick  response  in  England,  due  in  part 
to  its  enthusiastic  advocates,  J.  R.  Young  and  A.  De  Morgan.  It  is 
today  a  strong  rival  to  Newton's  method  in  both  England  and  the 
United  States.  It  is  widely  used,  though  less  than  Newton's  method, 
in  Germany,  Austria,  and  Italy.  In  Fiance,  Newton's  method  has 
held  undisputed  sway.  Horner's  method,  when  combined  with  the 
theorems  of  Budan  and  Fourier  and  of  Sturm,  is  effective,  secure,  and 
expeditious. 

5.     SUMMARY 

1.  In  the  13th  century  Ch'in  Chiu-shao  (1247),  Li  Yeh  (1259), 
and  Chu-chieh  (1299)  employed  a  method  of  approximation  vir- 
tually the  same  as  Horner's  method.  In  the  17th  century  it  was  in- 
troduced into  Japan. 

2.  Paolo  Rufnni  (1804)  invented  a  similar  method  in  Italy,  which 
was  soon  forgotten. 

3.  Horner  published  his  method  in  18 19,  and  it  soon  became 
widely  used  in  England  and  later  in  the  United  States,  and  to  a  less 
degree  in  Germany,  Austria,  and  Italy.  It  was  never  adopted  in 
France. 


X 

GENERAL  SUMMARY 

i.  The  earliest  known  method  of  finding  the  roots  of  numbers 
and  of  solving  formal  equations  is  that  of  Single  False  Position. 
This  method  was  used  by  the  Egyptians,  the  Greeks,  and  the  Hindus. 

2.  The  method  of  averages  was  probably  used  by  Archimedes; 
a  definite  exposition  of  it  is  given  by  Heron  of  Alexandria.  Its  appli- 
cability was  restricted;  but  it  was  often  incorporated  into  other 
methods,  as  in  Chuquet's  Rule  of  Mean  Numbers. 

3.  Methods  of  Double  False  Position  for  finding  the  roots  of 
numbers  seem  to  have  been  used  by  Archimides;  a  specific  process 
was  described  by  Heron  of  Alexandria.  In  algorithmic  work  double 
false  position  was  used  by  the  Arabs  and  Leonardo  of  Pisa.  By  the 
Renaissance  writers  these  processes  were  extended  to  the  solutions 
of  numerical  higher  equations.  Though  eclipsed  by  the  methods  of 
exhaustion  invented  by  Vieta,  Newton,  and  Horner,  methods  of 
double  false  position  have  never  become  obsolete.  Two  main  plans 
for  finding  mean  values  have  been  followed:  (1)  the  averaging  or 
analogous  devices,  as  illustrated  by  Chuquet  and  DeLagny;  (2)  in- 
terpolation by  proportional  parts  as  used  by  Leonardo  of  Pisa,  al- 
Karchi,  Biirgi,  and  Pitiscus. 

4.  Methods  of  exhaustion  have  been  used  the  most  extensively 
of  all  approximation  methods.  The  earliest  method  of  exhaustion 
recorded  is  that  used  by  Theon  of  Alexandria  for  irrational  square 
roots  by  reversing  the  binomial  formula.  By  a  special  adaptation  of 
this  method  Aryabhata,  Brahmagupta,  Sridhara,  and  Bhaskara 
worked  out  a  process  for  evolving  the  roots  of  numbers  digit  by  digit 
(Evolution).  This  was  adopted  by  the  Arabs  and  by  them  trans- 
mitted to  Christian  Europe  through  the  writings  of  Leonardo  of  Pisa 
and  of  Sacrobosco.  Oui  present  standardized  process  of  finding  roots 
of  numbers  is  largely  the  work  of  Peurbach  and  Cardan. 

Vieta  extended  this  method  to  numerical  higher  equations.  His 
method  was  popularized  and  made  more  practicable  by  Harriot, 
Oughtred,  and  Wallis,  and  remained  the  prevailing  method  until 
superseded  by  Newton's  method. 

5.  Newton's  method  of  exhaustion,  based  on  equations  whose 
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terms  constitute  an  infinite  series,  was  conceived  as  an  aid  to  his 
work  in  the  calculus.  The  place  of  this  method  as  a  constituent  part 
of  algebra  is  largely  due  to  its  modification  and  systematic  develop- 
ment by  Raphson.  Its  technique  was  improved  by  DeLagny, 
Halley,  Taylor,  and  Simpson  and  its  scientific  basis  was  clarified  and 
strengthened  by  Lagrange,  Fourier,  Budan,  and  Sturm.  No  other 
method  of  approximation  has  come  up  to  it  in  general  popularity. 

6.  A  third  method  of  exhaustion,  employing  continued  fractions, 
was  invented  by  Lagrange.  It  is  a  theoretically  perfect  method,  but 
too  laborious  for  practical  use.  It  has  retained  chiefly  an  academic 
interest. 

7.  Rolle's  Method  of  Cascades  and  the  method  of  recurring 
series  developed  by  Bernoulli  and  Euler  have  been  used  only  slightly. 
The  methods  invented  by  Collins,  Fontaine,  and  some  of  those  used 
by  DeLagny  were  never  taken  up  into  the  mathematical  activities 
of  the  world. 

8.  A  fourth  method  of  exhaustion  in  solving  numerical  higher  equa- 
tions is  both  the  oldest  and  the  youngest, — oldest,  because  it  was 
used  by  the  Chinese  in  the  thirteenth  century,  youngest,  because  it 
was  re-invented  and  given  to  the  Western  world  by  Ruffini  and 
Horner  in  the  twentieth  century.  It  is  a  process  of  evolution  of 
digits  resembling  that  of  Vieta;  however,  the  digits  are  evolved,  not 
by  the  inversion  of  a  binomial  formula,  but  by  the  transformation  of 
equations.  Horner's  method  has  made  progress  in  Italy,  Germany, 
and  Austria,  but  is  hardly  ever  used  in  France;  in  England  and  the 
United  States  it  is  a  close  iival  of  Newton's  method. 


APPENDIX 
Note  i 

The  solution  of  numerical  higher  equations  by  means  of  graphs  will  often 
give  as  near  an  approximation  as  is  desired.1  Where  greater  accuracy  is 
desired  a  combination  of  the  graphic  method  with  one  of  the  standard  methods 
of  exhaustion  or  of  double  false  position  is  generally  advantageous.  The 
application  of  the  calculus,  such  as  is  ordinarily  used  in  curve  tracing,  facili- 
tates the  work  of  approximation;  in  isolating  the  roots,  for  example,  one  can 
often  use  Rolle's  theorem  instead  of  the  more  elaborate  theorems  of  Fourier, 
Budan,  and  Sturm. 

Whether  tables  will  ever  be  used  in  the  solution  of  numerical  equations  is 
for  the  future  to  tell.  DeLagny  advocated  their  use  and  worked  out  a 
comprehensive  set  of  tables  similar  to  our  modern  integration  tables  in  his 
Analyse  Generale;  but,  as  his  voluminous  work  has  been  very  little  known,  his 
idea  has  never  been  tried  out. 

Note  2 

Teachers  of  algebra  often  fail  to  recognize  how  extensively  the  methods  of 
double  false  position  are  used  in  mathematical  thinking.  Notice  how  similar 
to  Burgi's  method  of  solving  equations  is  the  following  common  method  of 
determining  the  autumnal  equinox:  As  given  in  the  Almanac  of  1899,  the 
sun's  declination  at  approximate  noon  (Greenwich),  Sept.  22,  was  +  o°  18' 
8"  .  7  and  at  noon,  Sept.  23,  it  was  —  o°  5'  13"  .  9.  The  difference  in  declina- 
tion being  23'  22"  .  6,  the  equinox  is  found,  by  interpolation,  to  be  18'  18"  .  7 
4-  23'  22"  .  6  or  0.7762  of  a  day  after  noon,  Sept.  22.  Take  also  the  Law  of 
Charles  in  physics,  by  which  we  know  that  if  the  temperature  of  a  gas  be 

lowered  1°  C.  its  volume  will  be  decreased  by  —  of  its  volume  at  o*  C,  and  if 

273 

its  temperature  is  lowered  io°  C.  its  volume  will  decrease  by  —  of  its  volume 

273 

at  o°  C.  From  these  two  positions  the  volume  at  any  temperature,  say 
—  ioo°  C,  for  which  the  gas  does  not  become  a  liquid,  is  calculated.  The 
Fahrenheit-Centigrade  thermometer  problems,  the  interpolations  for  trigo- 
nometric functions,  and  many  problems  in  proportion  are  more  familiar  illus- 
trations. 

•  Consult  Arthur  Schultze,  Graphic  Algebra.  New  York,  1908. 
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The  formal  and  informal  interpolation  processes  used  in  astronomy,  physics, 
trigonometry,  and  other  computational  fields  outside  and  inside  the  class- 
room constitute  a  big  share  of  the  actual  calculation.  There  seems  to  be  no 
valid  reason,  then,  why  the  student  of  algebra  should  not  have  his  attention 
called  to  the  use  of  such  methods  in  the  solution  of  numerical  higher  equa- 
tions. While  less  refined  than  methods  of  exhaustion,  they  are  very  usable  for 
securing  results  correct  to  one  or  two  decimal  places.  An  interesting  text- 
book exposition  is  that  given  by  Professor  Isaac  Dalby  in  the  second  volume 
of  his  A  Course  of  Mathematics,  designed  for  the  use  of  the  Officers  and  Cadets  of 
the  Royal  Military  College  (London,  1806),  pp.  175-185.  For  good  illustrations 
of  the  use  of  these  methods  in  recent  texts,  see  Wells's  University  Algebra 
(Boston,  1884),  pp.  417-19,  and  F.  L.  Griffin's  Introduction  to  Mathematical 
Analysis  (New  York,  192 1),  pp.  37-38. 

To  find  the  mean  between  the  two  false  positions  we  generally  use  one  of 
the  following  three  ways:  (1)  finding  the  arithmetic  average,  or  using  some 
such  special  formula  as  Chuquet's  Rule  of  Mean  Numbers;  (2)  interpolating 
directly     by     proportional     parts;    (3)    using     the     interpolation     formula 

o»ei  —  axe2 

x  =  , where  a\  and  a2  are  the  approximate  values  and  e\  and  e2 

ei  -  e2 

are  the  corresponding  errors. 

Note  3 

The  solutions  of  numerical  higher  equations  first  became  a  mathematical 
need  in  the  work  of  geometry,  trigonometry,  and  astronomy.  However,  they 
enter  into  expert  practical  calculation  more  frequently  than  the  average 
student  realizes,  and  the  author  feels  convinced  that  the  insertion  of  such 
problems  would  vitalize  the  methods  of  approximation  given  in  our  college 
algebras.  The  following  representative  problems  from  the  fields  of  physics, 
mechanics,  economics,  and  engineering,  selected  from  the  more  progressive 
text-books,  will  serve  to  illustrate:2 

(1)  The  real  root  of  the  equation  x5  —  38.V  —  101  =  0  gives  the  diameter 
of  a  water  pipe  200  feet  long,  discharging  100  cubic  feet  per  second  under  a 
head  of  10  feet. 

Result:  diameter  is  2.92  ft. 

(2)  The  equation  -fi  —  192001*  +  21 1200  =  o  gives  the  speed  in  feet  per 
second  on  a  i-inch  manila  rope  transmitting  4  horse-power,  under  a  tension 
of  300  pounds  on  the  tight  side. 

:  Merriman  and  Woodward,  Higher  Mathematics,  p.  13;  Rietz  and  Crathorne,  Col- 
lege Algebra,  p.  140;  Karpinski,  Benedict,  and  Calhoun,  Unified  Mathematics,  p.  404; 
F.  L.  Griffin,  Introduction  to  Mathematical  Analysis,  pp.  37-42,  337-342.  More  prob- 
lems of  a  similar  nature  may  be  found  in  these  same  references. 
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Result:  speed  is  11.07  feet- 

(3)  From  the  American  Report  on  Wholesale  Prices,  Wages,  and  Trans- 
portation, for  1891,  the  median  wage  in  dollars  is  given  by  j^  of  a  value  of  x 

1  19  s  4^  107  S 

in  the  equation  2561  -  =  6972 \-  657  —  x  —  33  —  x2  H x3 x*. 

2  128  18  64  48  128 

Result:  wages  were  $1,536. 

(4)  In  Archimedes's  problem  of  the  section  of  a  sphere,  if  k  represents  the 
ratio  of  the  larger  to  the  smaller  segment,  the  distance  of  the  plane  section 
from  the  center  of  a  sphere  10  inches  in  diameter  is  given  by  the  equation 


300.V  +  2000 
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(5)  A  sphere  of  yellow  pine  1  foot  in  diameter  floating  in  water  sinks  to  a 
depth  x  given  by  the  equation  2x3  —  3X2  =  0.657. 

Result:  depth  is  0.606  foot. 

(6)  The  equation  d*  +  320J'2  —  340^  —  4290  =  o  gives  the  smallest  safe 
diameter  (in  inches)  for  a  certain  steel  shaft. 

Result:  the  diameter  is  4. 11  inches. 

(7)  The  equation  2.V3  —  9.Y2  -+-  24  =  o  gives  the  depth  of  immersion  in 

Q 

water  of  a  sphere  3  feet  in  diameter  and  with  a  specific  gravity  -. 

Result:  depth  is  2.38  feet. 

(8)  The  equation  x4  —  n8x2  +  i6ox  —  115  =  o  gives  the  length  of  the 
longest  rectangular  panel  1  foot  wide  which  can  be  fitted  diagonally  across  a 
door  4  feet  wide  and  10  feet  long. 

Result:  length  is  10.17  feet. 

(9)  The  "index  of  correlation"  between  the  eye-colors  of  a  certain  group  of 
people  and  of  their  great-grandparents  is  approximately  the  root  of  the  equa- 
tion 0.024.x4  -f  0.137X3  +  0.035*?  +  x  ~  °-225  —  o  which  lies  between  o 
and  1.    [C.  B.  Davenport] 

Result:  index  is  0.22. 

(10)  Using  Jupiter's  radius,  45090  miles,  as  the  unit  of  distance,  the  equa- 
tion x3  —  5xa  -f  6.27396.V  —  0.060385  =  o  gives  the  greatest  and  least  dis- 
tances of  Jupiter's  Fifth  Satellite  from  the  center  of  the  planet. 

Results:  distances  are  2.48  and  2.51  units. 

(11)  A  magnet  placed  with  its  ends  in  a  "magnetic  meridian"  will  neutralize 
the  earth's  magnetism  at  certain  points.    To  calculate  the  position  of  these 

20000X 

points  in  a  certain  case,  it  was  necessary  to  solve  the  equation =  02 

(x2  -  ioo)2 
Result:  x  =  47.82  and  0.09. 
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Note  4 

Bearing  upon  source  material  relating  to  the  solution  of  numerical  equa- 
tions, the  reader  may  be  interested  to  know  that  in  New  York  libraries  are 
found  certain  works  very  rare  in  this  country.  Among  these  may  be  men- 
tioned the  manuscript  copies  of  the  Lilawati,  the  Kholasat-al-Hisdb,  and  of 
Chuquet's  Tri party  in  the  private  collection  of  Professor  David  Eugene 
Smith.  In  the  Columbia  University  library  there  is  a  manuscript  on  algebra 
by  Scheubel,  in  his  own  handwriting,  prepared  for  the  press  but  never  pub- 
lished; also  a  copy  of  Vieta's  collected  works  published  by  van  Schooten  in 
Leyden  in  1646,  and  a  copy  of  J.  Hume's  exposition  of  Vieta's  algebra, 
published  in  Paris  in  1636.  The  latter  work  seems  to  have  been  unknown  to 
the  mathematical  world,  as  we  find  no  references  to  it  in  any  published  history 
of  mathematics.3  Hume  anticipated  the  "Newtonian  divisor"  by  an  analogous 
device,  but  his  work  was  evidently  unknown  to  his  contemporaries.  In  the 
New  York  Public  Library  are  found  three  comparatively  rare  books:  Heri- 
gones'  Cursus  Mathematicus  (Paris,  1642),  which  popularized  Vieta's  method 
in  France;  Michel  Rolle's  Demonstratio  d'une  Methode;  and  Raphson's 
Analysis  Aequationum  Universalis,  edition  of  1697,  bound  with  his  almost 
wholly  unknown  metaphysical  treatise,  De  Spatio  Reali. 

5  It  is  referred  to  in  Professor  D.  S.  Smith's  forthcoming  history  of  mathematics. 
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